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XXXVITI. Isotopic Spin Selection Rules 
VII : Breakdown of the Rules and the Situation in %O 


By D. H. Witxinson 
Cavendish Laboratory, Cambridget 


[Received January 26, 1956] 


ABSTRACT 


Factors governing the breakdown of the isotopic spin selection rules 
are discussed. Two energy ranges exist within which we expect the 
rules to be obeyed ; at low excitation individual excited states are well 
separated in energy and have good isotopic spin purity so that the rules 
are obeyed in reactions passing through such states only; at high 
excitation individual states no longer have a well-defined isotopic spin 
but the rules are obeyed because many overlapping levels are simul- 
taneously excited in a total state of initially well-defined isotopic spin 
which breaks up before the Coulomb forces have had time to cause much 
isotopic spin mixing. In between these regions where the rules are 
obeyed (although for very different reasons) severe breakdown is to be 
expected. The approximate extent of these reasons is discussed. The 
experimental situation as regards the operation of the rules and other 
consequences of charge independence is examined in !%O where we 
have knowledge relating to the effective isotopic spin purity of ten 
energy levels or energy ranges between 6 and 30 Mev excitation. Good 
correspondence with theoretical expectation is found in the above- 
described general trend of effective impurity versus excitation. It is 
concluded that the situation is adequately described by the Coulomb 
perturbation of an otherwise charge independent (symmetric) Hamiltonian : 


$1. INTRODUCTION 


Durine the past few years a body of data has built up concerning the 
validity of the concept of isotopic spin as applied to the classification of 
states of light nuclei. Attention has focused on the one hand on the 
comparison in mass between members of isotopic spin multiplets (see 
paper VIII of this series) and on the other on the operation of the selection 
rules associated with the isotopic spin quantum number. The first type 
of study is of chief interest in attempting to get a measure of the charge- 
independence of nuclear forces, in particular the relationship between nn. 
and n—p forces. The effect of the Coulomb perturbation in disturbing 
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the results to be expected under complete specifically nuclear charge 
independence is here a second order one: the first order effect is the 
mixing by the Coulomb forces of the isotopic spins of states of different 
primitive isotopic spin but of the same ordinary spin and parity; this 
is followed by the ordinary configuration mixing of the components of 
the same isotopic spin in the various sets of states of different primitive 
isotopic spin with attendant displacement in energy in the usual way. 
We may therefore hope that the Coulomb perturbation will not interfere 
drastically with investigations of charge-independence through the 
comparison of multiplets. On the contrary in the study of the selection 
rules we see the relaxation of the rules as a first order consequence of 
the Coulomb perturbation which is consequently much easier to pick out. 

We have in this series of papers been chiefly concerned with determining 
the degree to which the selection rules hold and with attempting to assess 
in each case whether or not the result was consistent with the expected 
effects of the perturbation of a completely charge-independent (symmetric) 
Hamiltonian by Coulomb forces alone. This is a difficult procedure 
because it is only for especially simple and usually experimentally 
inaccessible states that one can compute the isotopic spin mixing with 
any degree of confidence. The matrix element of the Coulomb perturba- 
tion between states of different primitive isotopic spin varies by more than 
a factor of ten from case to case so the theoretical relaxation of the 
selection rules which is proportional to the square of the mixing is often 
unsure to two orders of magnitude. However sufficient data are now 
becoming available in a few cases for us to begin to survey the situation 
as regards the relaxation of the rules—the isotopic spin impurity—over 
a wide range of excitation in a single nucleus and to see whether the 
general trend of impurity versus excitation that we observe is consistent 
with what we expect from the Coulomb perturbation. 

In this paper we attempt such an examination of 10. Before we do 
this we discuss in general terms the way in which we might expect the 
breakdown of the selection rules to depend on energy of excitation. 


§ 2. Errective Isotropic Spin Ivpurrry 


In discussing the degree to which we expect the isotopic spin 
selection rules to be obeyed we must distinguish two extreme classes 
of experimentation : 

(i) Experiments in which all the states that are involved from the 
isotopic spin point of view, including those of the intermediate compound 
system, are well-isolated and discrete. 

(ii) Experiments in which the states of the intermediate compound 
system whose isotopic spin is relevant to the selection rules are in the 
continuum region of strongly overlapping levels. 

And of course we shall in practice be interested in all intermediate cases. 
We first discuss these extreme classes : 
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(i) In this class are experiments such as those on the forbidden emission 
of 1 radiation with which we have chiefly dealt in this series. We deal 
only with well-isolated levels and so are concerned only with essentially 
static descriptions of those states such as a unique spin and parity and a 
time-independent isotopic spin make-up. In such cases in order to assess 
the severity of the expected breakdown of the rules we ask for «,(T7"), 
the relative amplitude of the state of isotopic spin 7’ that the Coulomb 
forces mix with our state of primitive isotopic spin 7'.+ If the process 
which we study is forbidden for isotopic spin 7 and allowed for 7’ we 
then expect it to proceed at a rate of order «,2(7T’) of that which we 
should have anticipated but for the isotopic spin inhibition (more properly 
the breakdown will depend on the appropriate combination of the «’s 
of the initial and final states between which the forbidden transition is 
made). Now «,(7’)=H%}7/4E pp where 4Ep, is the separation in 
energy of the state of 7 and those of 7’ of the same J,7 whose 
mixing by H*° we consider (we sum over all states 7’). Radicati 
(1953) and MacDonald (1954, 1955) have estimated H° in several cases 
in light nuclei. The treatment of MacDonald is rather the more detailed 
and we now make our discussion chiefly in the light of his results. 
a,(1’) may be very roughly broken down into two parts, one associated 
with the isotopic spin impurity arising from excitation of the core and 
the other from the mixing of the valence nucleons either with other 
configurations or with other states of their own configuration or mixed 
configuration. In the first part H° is large but so is 4H and «? due to 
core excitation of ls- or 1p-shell nuclei for example is about 2 to 4 107°. 
In the second part, when the contaminated state belongs to a pure 
configuration, the only case that MacDonald considers, H° is rather 
small as a rule, typically 0-1 Mev or less. Larger values up to nearly 
1 ey are to be found but only between states in which one of the particles 
changes its total quantum number without changing its orbital quantum 
number and such states are expected to be separated by 10-30 Mev so 
the large H° does not bring a large «*. For the more usual cases we may 
take H°=0-1 mey as a fairly safe upper limit although of course here 4H 
may perhaps be very small. We should also consider contamination 
arising in mixed configurations of valence nucleons. It seems that here 
we may expect both H° to be large and AEH perhaps to be small on 
occasion. It is well known (e.g. Lane 1954) that the interaction of an 
s- or d-wave nucleon with p-shell nuclei tends to be weak, giving rise to 
states whose parentage description in respect of that nucleon is rather 
simple. When, say, an s-wave nucleon interacts with the p-shell nucleus 


+ The £1 selection rule may also be violated by virtue of higher-order terms 
in the #1 Hamiltonian (Gell-Mann and Telegdi 1953) or of the inadequacy of 
the assumption that the charge resides on the protons—mesonic effects 
(Morpurgo 1954). Both these effects are small and are negligible for the 
following discussions. 
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(ground or excited state) J,z, 7’ we may expect them to form the closely- 
related states J-+4. 7, 74 with not very great separation which are able 
strongly to contaminate one another. 

So for the experiments of class (i) we should not expect to find the 
isotopic spin selection rules affording factors of discouragement greater 
than 500 or so since we always have an impurity of order 2 to 4x 107% 
from core excitation. When states of the same J, 7 and different 7’ 
approach within 500 kev or so of each other we must anticipate that 
valence nucleon contamination with H°~0-1 Mev may have taken over. 
Indeed at considerably greater spacing contamination of the mixed 
configuration type last discussed may have become considerable; by 
the time 14H=—1 mev has been reached the mixing may be almost complete 
and isotopic spin may have largely lost its meaning as a label for this 
state. This appears to have happened for 1°B at 6-88 Mev (Wilkinson 
and Clegg 1956). It will probably start happening for light (4 <20) 
4n+2-type nuclei at an excitation of 6 to 10 Mev, for light 4n-type nuclei 
at 14 to 22 Mev and at 11 to 17 Mey for light nuclei of odd mass. 

(ii) In this second class are experiments of the (d, d’), (d, «) and other 
types in which, although the initial and final states are well isolated and 
so receive descriptions of the time-independent type discussed above, 
the intermediate state may be at an excitation of perhaps 30 or 40 Mev 
in the continuum region of strongly overlapping levels. It was in such 
types of experiments as this that the applicability of the isotopic spin 
rules (straightforward conservation between initial and final states) was 
first established. However, at the high excitations of the intermediate 
state through which the reaction must pass and through which isotopic 
spin must be conserved, the considerations of (i) above would lead us 
to suppose that the mixing of isotopic spins of individual states should 
be so complete that isotopic spin should no longer have any significance. 
Why then is it that the simple conservation rules seem to work under 
these conditions ? The answer is that we are no longer dealing, as in (i), 
with single states in the intermediate compound system but with the 
simultaneous excitation of a great number of overlapping states of 
comparable amplitude. The correct picture is therefore no longer a 
time-independent one but one that involves the changing phase relation- 
ships between the many excited states. In particular the total state of 
the compound system as originally formed will approximate itself as 
closely as possible to the isotopic spin of the initial system. This is 
possible, despite the fact that the time-independent description of each 
individual state contains a mixture of several isotopic spins in comparable 
proportions, because we simultaneously excite many such states and a 
suitable choice’ of their phase relationships can result in an initially 


eee 


} Although it seems reasonable that large values of H¢ should be generated 
in this way they have not yet been found by direct calculation (J. P. Elliott— 
private communication), 
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more or less well-defined isotopic spin for the total state. The Coulomb 
forces perturb this total state and cause the growth of the other isotopic 
spins with a characteristic time of order 4/H°. If they had an indefinite 
time in which to act they would cause complete mixing and the system 
would forget the isotopic spin of the initial system from which it was 
made up. The final break-up would then lead to states violating the 
simple isotopic spin conservation rules as readily as to those obeying 
the rule. However, an indefinite time is not available and the underlying 
individual states have a finite width I which may be several Mev at the 
excitations we now discuss, so the break-up occurs with a characteristic 
time of order h/I’. If now I'S H° the system breaks up before the Coulomb 
forces have had time to cause forgetting of the isotopic spin of the initial 
system and the conservation rules are obeyed. The degree to which they 
are obeyed will now be determined by H°/I° rather than by H*/AE.+ 
So as we tend to higher and higher excitation the selection rules should 
be better and better obeyed as I” becomes bigger and bigger and more 
and more overlapping states are involved. A similar argument applies 
to reactions whose mechanisms are referred to as ‘direct’: again it is 
the short time scale that enables the conservation to hold. This class (ii) 
situation is characterized by a smooth excitation function with no clear 
remaining resonance structure. It will probably have set in by about 
14-18 mev in light 4n+ 2-nuclei, by 22-30 Mev in light 4n-nuclei and by 
15-25 Mev in light nuclei of odd mass. 

Between these two extremes, at intermediate excitations, we shall find 
situations where levels are dense and overlapping but where a sufficiently 
large number is not yet involved for the considerations of class (ii) to 
apply. Here we should expect to find severe breakdown of the selection 
rules. This intermediate situation will be characterized by an excitation 
function which although it is not one of well-separated resonances yet 
retains some bumps telling of traces of single level structure ; the total 
state is one of many overlapping levels but not of such a mixture that 
a few do not retain some individual importance ; we do not yet excite 
so great a mass that. we can arrange for the intermediate compound 
state to start off closely approximating the isotopic spin conditions of 
the initial system. 

At low excitation (class (i)), the operation of the selection rules is 
characterized by «2 the isotopic spin impurity of the individual states 
concerned. At high excitation (class (ii)) we have seen that the operation 
of selection rules is no longer determined by the properties of individual 
states, all of which are very impure, but rather by the properties of the 
assembly of overlapping states and by the way in which it is excited. 
We can still define an effective isotopic spin impurity «,_ as the reciprocal 
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; This is a rather picturesque way of describing the situation. A proper 
treatment would not split up the process into ‘ contamination ’ and ‘ decay 
but treat them together—with the same result. 
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of the factor by which the forbidden transition is inhibited remembering 
that it has nothing to do with a? of the individual states at that excitation 
but that it depends on many states and also on the way in which the 
assembly is excited. 

In general terms then we expect the operation of the selection rules 
to get worse with increasing excitation as «? increases, to pass through 
a region of more or less complete breakdown with values of «* ranging 
up to unity and finally to get better again as a,” takes over from «* 
in the continuum region. 


§ 3. SravTeEs oF 160 


We have no direct information relating to the isotopic spin purity of 
the ground state of 160 and for purposes of discussion of transitions 
involving it we adopt the theoretical figure of MacDonald, namely 
o92(L)~4 x 10-8, 

We now take the states in turn about which we have information 
relating to the purity of their isotopic spin and follow them with a 
discussion of the regions for which we can say something about o,°. 


(i) 6-14 and 6-91 Mev 

The observed mode of decay of the 6-91 Mev 2+7'=0 state? is by 
#2 radiation to the ground state rather than by #1 radiation to the 
6-14 Mev 3—7'=0 state. The systematics of gamma-ray transitions in 
light nuclei (Wilkinson 1953 a, 1955 a) suggest | W |?—=0-032 in Weisskopf 
units as the most probable a priori strength to expect for #1 transitions 
that do not violate the isotopic spin selection rule.t If we take | M |?=1 
for the competing #2 transition§ we should then have expected the 
relative strengths of the #1 and £2 transitions to have been 0-13: 1 
but for the violation of the isotopic spin selection rules. In (1) of this 
series Wilkinson and Jones (1953) reported that in fact this ratio is 


+ For the general assignments see Ajzenberg and Lauritsen (1955). 

i It is observed that 85°, of all allowed #1 transitions in light nuclei have 
speeds within a factor of 7 of | M |? =0-032. We shall use this figure when 
we wish to guess uninhibited transition strengths. Its reliability is then 
presumably rather better than an order of magnitude. 

§ There is little to guide us in guessing the strengths of #2 transitions 
Several examples of | |? >1 are known in light nuclei but so are some 
considerably lower values. It is unlikely that this particular #2 transition 
is very considerably enhanced since a very strong #2 transition, representing 
most of the available collective strength from the ground state, is known to 
the 14-7 Mev state (Wilkinson 1955 b). On the other hand Devons et al. (1955) 
have shown that the #2 lifetime is less than 1-2 < 10-14 seconds or | M |? >0-5. 
An upper limit is provided by the appropriate sum rule relating to E2 transitions 
that do not change the isotopic spin (Gell-Mann and Telegdi 1953). This 
limit is M |?~10 for this transition but is unlikely to be approached for the 
reasons just given. It therefore seems that | |? =1 is probably not too 
bad a guess, 
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<5x10: 1. Since then the 770 kev #1 transition has been sought 
by an independent (coincidence) method by the writer and Dr. B. J. 
Toppel at Brookhaven National Laboratory. This work confirms the 
earlier result and shows that the ratio is <2x10-3: 1. We therefore 
conclude, within the limits given by the uncertainty of our guesses about 
uninhibited widths, that the isotopic spin rule is imposing an inhibition 
factor of at least 65 on this transition and so «,2(1) for both of these 
levels is <0-015. ) 
(ii) 7-12 Mev 

Devons, Manning and Bunbury (1955) report that the lifetime of 
the 1—7'=0 level for H, transition to the ground state is less than 
8x10-'° seconds. This means | M |?>2-3x10-4. Since the ground 
state 7’=1 impurity alone would produce | M |?~10-4 we conclude 
that probably for the excited state «2(1)>4 x 10-8 (neglecting interference 
effects). 

(ili) 8-87 Mev 

This 2—7'=0 state was observed by Toppel, Wilkinson and Alburger 
(1955) to decay to the 2+7'=0 state at 6-91 Mev with about one-third 
the intensity of its decay to the 1—7'=0 state at 7-12 Mev. This would 
appear to suggest the high 7’=1 impurity for the 2— state of «92(1)~0-03+ 
or a little less if allowance has to be made for the impurity of the 
6-91 Mev state—see (1) above. However, it was also observed indirectly 
that the M1 competing transition is very probably rather slow and it 
was concluded that in fact this forbidden #1 transition is more probably 
inhibited by a factor of about 900 (see Toppel et al. 1955 for the detailed 
argument). So we say 10-§<«,?(1)<0-03 with the lower figure favoured. 


(iv) 12-95 Mev 

This 2— state is probably the first T7=1 state of 1°O (see VIII of 
this series). However it decays by alpha-particle emission to the 
2+ 4-43 mev 7'=0 state of #2C with a reduced width of 4-2% (alter- 
natively 11-:2°%) of the single-particle value.{ Although the single- 
particle value is not an absolute upper limit it has never been surpassed 
with certainty. It is of course practically the same as the Wigner sum-rule 
limit. An alpha-particle reduced width of 4% is quite a large one even 
for an allowed transition; we can therefore say «,7(0)>0-04 with the 
likelihood that it may be considerably larger. The 2— 7=0 state at 
12-51 Mev resembles this present one in having a very large reduced 
width for d-wave proton emission to the ground state of 1°N (75% of the 
ee 

+ In making this estimate we use the result from the systematics that the 
most probable | |? value for M1 transition in light nuclei is about 0-15 
(an uncertainty of roughly an order of magnitude either way covering most 


of the transitions). 
tWe use a=1-45(A,¥3+4,"3)x 10-4 cm, a recipe found suitable for 


the analysis of alpha-particle scattering by C (Hill 1953). 
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single-particle value against 45°, for the present state). It is therefore 
possible that these two states are the 2— 7’=O and 1 linkages of a 
1d5,, proton to the ground state of 1°N as unique parent. They would 
then contaminate each other rather effectively in the manner discussed 
in § 2 for mixed configurations. In this case our alpha-particle emission 
from the 7’=1 state is due to the admixture of the 70 character from 
the other, whose reduced alpha-particle width is 80°, of the single-particle 
value. If this interpretation is correct we obtain the value «,?(0)~0-05 
for the present state ; in case it is not correct we must continue to use 
0-04 as a lower limit. 
(v) 13-05 Mev 

This 1— state behaves like 7’=1 in that its #1 transition to the ground 
state shows | M |?~0-1, a large value. But it also decays by alpha- 
particle emission to the 0-+ ground state and 2+ 4-43 Mev state of 2C, 
both 7'=0, with reduced widths of 2-8 and 5-5°%, (alternatively 5-6 and 
11-0°%) of the single-particle value respectively. So we say «,7(0)>0-05 
with the probability that this is a low figure. The #1 radiative width 
for the ground state transition is 150 ev (Schardt, Fowler and Lauritsen 
1952). The large s-wave reduced proton width to the ground state of 1©N 
agrees well with that theoretically expected if the ground state is the 
unique parent. In this case we can with some confidence compute the 
expected radiative width in either extreme coupling scheme under the 
assumption of charge imdependence (a pure 7’=1 state). We find 
theoretical widths of 150ev in LS coupling and 50 ev in jj coupling. 
The good agreement with experiment found in LS coupling which is 
expected to be the better approximation suggests that in fact the 7’=0 
contamination cannot be complete. «,?(0)=0-5+ would reduce the 
theoretical radiative width to about 100 ev or less which begins to 
conflict with experiment. We therefore say 0-06<«,2(0)<0-5. The 
relevant spin and parity assignments have been made earlier (see 
Wilkinson 1953 b); since then the work of Kraus (1954) has confirmed 
them. 

(vi) 16-19 Mev 

In this region resonance structure has been observed and measured 
in the reaction 1®O(y,n)%O by Katz, Haslam, Horsley, Cameron and 
Montalbetti (1954) and by B. M. Spicer and A. Penfold (unpublished work 
at the University of Ilinois—see also Ajzenberg and Lauritsen 1955). 
Resonance structure in the same region has been observed in 16O(y, p)15N 
by Johansson and Forkman (1955). This region of excitation is low 
enough to make it very likely that only the ground state of 15N is involved 
in the (y, p) reaction and only the ground state of 1°O is energetically 
accessible in the (y,n) reaction. If the isotopic spin is well defined we 


} Such large values of «* are of course no longer relevant to any perturbation 
expansion but merely indicate relative intensities for the components of 
different isotopic spin. 
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then expect the neutron and proton yields to bear calculable ratios to 
each other; these ratios do not depend critically on the (unknown) 
multipolarity and type of transition. In general the (y, p) resonances 
are found at the same positions as the (y,n) and there is also a quite 
good general correspondence between the experimental relative intensities 
and those expected under charge independence. This leads to the 
conclusion that in this region isotopic spin is a fairly good quantum 
number. If we take «? as a measure of the fractional deviation from 
the expected intensity ratios we say «?<0-3. This figure is obtained by 
clumping the observations in units of 1 Mev. If we compare individual 
transitions some deviations from expectation are apparent that would 
lead us to say «?>0-1 or more for some levels but the statistical accuracy 
is very low and there are discrepancies between the various measurements 
of the level strengths in the (y, n) process used for the comparison so 
we do not feel safe in concluding that even this small amount of impurity 
is definitely established. 


(vil) 20-23 Mev 

This is the region of the ‘ giant resonance’ seen in 1%O(y, n)!5O at 
22-2 Mev and 1%O(y, p)4°N at 22-3 Mev. The structure of these resonances 
has been investigated by Katz et al. and by Spicer and Penfold in the 
first reaction, and by Stephens, Mann, Patton and Winhold (1955) in 
the second. There is good general agreement as to position and gross 
features of the resonance seen in the two reactions and this we should 
expect whether the isotopic spin were well-defined or not. There appear, 
however, to be some features not consistent with well-defined 71 (the 
resonance is #1). For example, the ratio of the integrated cross sections 
in the ranges 21-22 Mev and 20-21 Mev is about 1:1 for (y,n) and 
25:1 for (y,p). (Since the binding energy of a neutron in 1%O is 
16-0 Mev and the first excited state of 15O is at 5-27 Mev we can ignore 
transitions to all but the ground state of 1°O in the giant resonance 
region.) The figures given above for (y, p) refer to emission to the 
ground state of 4*N only. Again the integrated cross section to the 
peak of the resonance in (y,n) is about 20 mMev-mb whereas in (y, p) 
it is only about 14 Mev-mb (again to the ground state only). This is 
despite the fact that the proton energies are well above the Coulomb 
barrier and nearly 4 Mev more energy is available to the protons than 
to the neutrons ; this would lead us to expect roughly 30° more protons 
than neutrons. On the basis of these observations we should suppose 
that the requirements of charge-independence were no longer being 
fulfilled even though no actual selection rule was being broken. This 
degree of difference of behaviour for neutrons and protons would seem 
to imply at least «,?(0)>0-1. 

We must also consider, however, the behaviour of the (y, ~) and allied 
reactions in this region. These have been analysed by Gell-Mann and 
Telegdi (1953) and by Hsiao and Telegdi (1953) as suggesting that the 
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E1 cross section for '6O(y, «)!2C is suppressed in the region of the giant 
resonance because the resonance state is of 7=1 while a T=1 state 
in 2C is not yet energetically available (allowance being made for the 
Coulomb barrier). 7'=1 states in 2C become available above L,,~25 Mev 
Above this energy the cross section for (y, «) remains at more or less the 
same level as in the giant resonance region and transitions are made 
chiefly to the 7=1 region in #C. It has been supposed that this indicates 
a considerable degree of isotopic spin purity in the giant resonance region 
around 22 Mev since if the total absorption cross section were dropping 
rapidly from the resonance peak as we move to higher energies the 
constancy of the (y,«) cross section would indicate that mechanisms 
suppressed in the giant resonance region (H1 absorptions) were coming 
into play above the effective threshold for transition to 7’=1 states 
in 20, However, we do not believe that this inference can necessarily 
be made since in fact the total absorption cross section at H,—25 Mev 
seems to be fully a third of its peak value according to Stephens et al. 
(compare the ‘long tails * observed in ?C(y, n)!4C and 1O(y, n)?®O by 
Sagane (1951) and in more detail by Barber, George and Reagan (1955) 
in 2C(y,n)"4C). The cross section for 1*O(y, 4x) has also fallen somewhat 
from 22-6 Mev (where it too has a resonance) to 25 Mev (Goward and 
Wilkins 1952) and so it has not changed its form dramatically relative 
to the total absorption cross section. Indeed the very rapid rise in the 
16O(y, 4) cross section above 20 Mev which coincides with the rapid 
rise in the total cross section and the coincidence of its strongest 
resonance with the giant resonance suggest that in fact #1 absorption 
is taking place for the (y, ~) processes in the giant resonance region and 
that H2, M1 and other mechanisms are relatively unimportant there. 
We therefore believe that the evidence from the (y, «) reactions is quite 
consistent with our conclusion that in the giant resonance region there 
is considerable isotopic spin impurity. On the other hand if we go to 
somewhat higher energies (H,,>30 Mev) the (y, «) processes have increased 
by a factor of 4 or 5 in strength relative to the (y, n) and this suggests 
that they are here enjoying some freedom denied them at the lower 
excitation where therefore «?=1 is presumably not attained. It seems 
then that all evidence is consistent with a value of «,2(0) of about 0-1—0-5 
in the giant resonance region. 


(vili) 25-35 Mev 

Above 25 Mev we continue the argument based on the !O(y, «)!2C 
reaction. Hsiao and Telegdi find that here the mechanism is one that 
appears to imply a fair degree of effective purity to the isotopic spin 
states. This they see because transitions to a (7’'=0) state in 2C at 
about 13 Mev which were strongly favoured at lower energies give place 
to transitions chiefly (by a factor of 5 or more) to the 7’=1 region of 12C 
around 16 Mev. This then suggests «?(0).~<0-1 in this region. Consistent 
with this improvement in effective isotopic spin purity that has taken 
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place above the giant resonance region is our observation in (vi) above 
on the trends of relative (y, ~) and total absorption cross sections. 


(ix) 26-8 Mev 

This is the energy of the compound system when the first: 77=1 state 
of *4N at 2-31 Mev fails to show up in 14N(d, d’)!4N at Hqg—6-98 Mev 
(Bockelman, Browne, Buechner and Sperduto 1953). The 7'—1 state is 
excited at least 20 times less strongly than the 70 state at 3-95 mev 
so we say for the highly excited 1°O state x 2(1).—~<0-05. We must note 
in caution, however, that although this 7'=1 state is strongly excited in 
14N(p, p’)4N at Hp=6-92 mev—half as strongly as the 3-95 mev state 
(Bockelman et al.)—it is only weakly excited at H)—9-6 Mev (Freemantle, 
Gibson, Prowse and Rotblat 1953). ; 


§ 4. Discussion 


We will first of all decide in the light of § 2 what values of «? or a4? 
we should expect in this range of excitation in 1°O. Before levels begin 
to overlap very strongly we have three types of contamination to 
consider : 


(i) that due to core impurity ; 
(ii) that due to mixing of valence nucleon states when the contaminated 
state is of a pure configuration ; 
(iii) as (ii) when the contaminated state is not of a pure configuration. 


Type (i) contamination persists throughout and may be presumed to 
lie between 10-3 and 10-2. Type (11) contamination we also find every- 
where and it may be dominant when we are remote from the rather 
special levels which may have very heavy contamination of type (iii). 
But in the present case there is no type (ii) contamination proper because 
16Q) is a closed-shell nucleus. However we may anticipate that type (iii) 
contamination will also show two forms : the first belongs to the special 
states of unique parentage discussed in § 2 where H° may be large ; the 
second belongs to the everyday states with complicated parentage 
relationships where smaller values of H° similar to those of type (ii) 
contamination will be found. 

The first 7’=1 state in 180 is probably that at 12-95 Mev so until we 
approach within a few times H® of this energy the impurity will remain 
at the type (i) level of 10-3 to 10-2. We must then expect that type (i) 
configuration may make itself felt above about 9 Mev since it is long 
range. For example the 2— 7'=0 state at 8-87 Mev may be contaminated 
by the 2— first 7'=1 state at 12-95 mev. In fact we know that these 
two states are quite closely related because 1°N shows an allowed 
beta-transition of log ft=4:4 to the 8-87 Mev state (Toppel ef al. 1955) 
and so H* may well be quite large. So above an energy of 9 Mev or so 
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x2 may begin to rise rather sharply and reach a value of order unity by 
about 12 Mev (complete isotopic spin mixing). On the low sidé we 
anticipate a rise due to the ‘ type (ii) ’ aspect of type (ili) contamination. 
At 13 mev states of the same spin and parity are separated by about 
Imev. Taking H‘~0-05 Mev as a reasonable lowish figure for this type 
of mixing we see that we may anticipate «?~3 or 4x10 at 13 Mev 
due to this mixing and core impurity combined. As the excitation 
increases above the onset of 7’=1 states the level spacing diminishes 
until at 22 Mev it is down to about 200 kev or less with a corresponding 
lower limit to likely «? values of about 0-05. 

At about 22 mev, although some levels are still sharp, broad levels 
are common and we approach the region of complete overlap where «? 
is replaced by «2 as we discussed in §2. By the time an energy of 
30 Mev is reached it seems likely that good overlap will be taking place 
with [~1 Mey or more. The larger values for H° are no longer relevant: 
since we deal with the properties of the large numbers of levels into which 
single-particle motions will now be broken. On the other hand the 
smallest values of H°® need no longer be considered since their effects will 
be swamped by the stronger perturbations. H°~0-1 Mev now seems a 
reasonable value giving «,_?~0-01. The upper and lower limits will 
converge as the excitation energy increases and we tend to a completely 
statistical description (which will, however, depend on the mechanism of 
excitation). 

We can now compare these theoretical expectations with the results of 
experiment summarized in §3. On the figure we display as the two full 
lines the theoretical rough upper and lower limits to the effective isotopic 
spin impurity that we have just derived following § 2. We also show the 
experimental results. The meanings of the various ways of entering the 
results will be clear from the discussion in § 3. The three lowest energy 
entries are ‘limits’ (within the uncertainties involved in guessing the 
uninhibited radiative widths). The next entry, at 8-87 Mey, has an 
upper limit (in the same sense) but we have evidence that the true 
impurity is much lower so we enter the point at the low value and 
indicate by the downwards pointing arrow that we attach more weight 
to the location of the point than to the upper limit. The next entry 
has a fairly strict though not inviolable lower limit given by the single- 
particle value and which we indicate, but we enter the point at the value 
suggested by the identification of the contaminating state as that at 
12-51 Mev, and indicate by the broken arrow that we think this identi- 
fication a likely one (i.e. that we probably have here an actual value for «2 
rather than a limit only). For the next entry we have a fairly firm lower 
limit at which we place the point, indicating by the upwards pointing 
arrow that the impurity could run as high as shown but probably not 
all the way to #*=1. The next entry is at 16-19 Mev and is an upper 
limit, a lower limit being shown dashed because we do not regard it as 
well established. The next entry has both upper and lower limits that 
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seem equally well established so we place the point midway between 
them. The entries at 27 Mev and 25-35 mev are again upper limits only. 


IMPURITY 


O1 
Ls) 


ISOTOPIC SPIN 


| 
) 


lO 


BREE GIVE 


OF lO ; 20 30 40 
EXCITATION ENERGY MEV 


The theoretical and experimental isotopic spin impurities in 1O as a function 
of excitation. The two full lines are the upper and lower limits 
theoretically expected for the impurity (defined as ao? or ag? as is 
appropriate). The experimental points have the significance discussed 
in the text. Those around 7 Mev and 13 Mev have been shifted sideways 
a little to avoid overlapping. 


Although the theoretical estimates and the extractions of the impurity 
from the experimental data are both admittedly rough it is clear from 
the figure that there is a very good agreement in general terms between 
the two. In particular we see in the experimental situation the expected 
high purity at low excitation followed by very substantial impurity at 
medium excitation and then by a return to a good effective purity at 
high excitation (for very different reasons) as a," takes over from «?. 


§ 5. CONCLUSION 
We conclude that the isotopic spin properties of 1°O over a wide range 
of energy show considerable changes with excitation but are in good 
agreement with theoretical expectation based on the Coulomb perturbation 
of an otherwise charge-independent (symmetric) Hamiltonian. 


392 Isotopic Spin Selection Rules: VII 


ACKNOWLEDGMENTS 


I should like to thank Dr. B. J. Toppel for his participation in the 
measurement described in § 3 (i). The discussion in § 2 of the operation 
of the selection rules in class (ii) experiments has benefited from 
conversations with many people particularly Dr. A. M. Lane, Dr. J. 
Weneser, Professor D. C. Peaslee and Professor R. F. Christy. To all 
these gentlemen I am pleasurably indebted. 


REFERENCES 


AJZENBERG, F., and Lauritsen, T., 1955, Rev. Mod. Phys., 27, 77. 

BarBeErR, W. C., GeorckE, W. D., and Ruaaan, D. D., 1955, Phys. Rev., 98, 73. 

BockELMan, C. K., Browne, C. P., BuEcHNER, W. W., and SPERDUTO, A., 
1953, Phys. Rev., 92, 665. 

Devons, S., Mannine, G., and Bunsury, D. Sr. P., 1955, Proc. Phys. Soc. A, 
68, IS. 

FREEMANTLE, R. G., Grsson, W. M., Prowsz, D. J., and Rorsuat, J., 1953, 
Phys. Rev., 92, 1268. 

GELL-Mann, M., and Trevzent, V. L., 1953, Phys. Rev., 91, 169. 

Gowakp, F. K., and Wixxtns, J. J., 1952, Proc. Phys. Soc. A, 65, 671. 

Hiri, R. W., 1953, Phys. Rev., 90, 845. 

Hsrao, C. A., and Teieenpt, V. L., 1953, Phys. Rev., 90, 494. 

Jonansson, 8. A. E., and Forxkman, B., 1955, Phys. Rev., 99, 1031. 

Katz, L., Hastam, R. N. H., Horstzy, R. J., Cameron, A. G. W., and 
MonTALBetti, R., 1954, Phys. Rev., 95, 464. 

Kraus, A. A., 1954, Phys. Rev., 94, 975. 

Lane, A. M., 1954, A.H.R.E., Document T/R 1289. 

MacDona.p, W. M., 1954, Thesis, Princeton ; 1955, Phys. Rev., 100, 51. 

Moreurao, G., 1954, Nuovo Cim., 12, 60. 

Raopticatt, L. A., 1953, Proc. Phys. Soc. A, 66, 139; 1954, Ibid. A, 67, 39. 

Sacang, R., 1951, Phys. Rev., 84, 587. 

ScuaRpt, A., Fow.er, W. A., and Lauritsen, C. C., 1952, Phys. Rev., 86, 527. 

STEPHENS, W. E., Mann, A. K., Patron, B. J., and WrnnHoLp, E. J., 1955, 
Phys. Rev., 98, 839. 

Topre., B. J., Witxinson, D. H., and ALBuRGsR, D..E., 1955, Phys. Rev., 
99, 632 (and in press). 

Wixmnson, D. H., 1953 a, Phil. Mag., 44, 450; 1953 b, Phys. Rev., 90, 721; 
1955 a, Phil. Mag. (in press), A.8.C.L. Chalk River Document PD-260 ; 
1955 b, Phys. Rev., 99, 1347. 

Wickinson, D. H., and Ciece, A. B., 1956, Phil. Mag. (in press). 

Witxinson, D. H., and Jonzs, G. A., 1953, Phil. Mag., 44, 542. 


fac Os ee! 


XXXVIII. Diffusion of the Chloride Ion in NaCl 


By D. Patrserson,} G. 8S. Rosz,+ and J. A. Morrison 
Division of Pure Chemistry, National Research Laboratories, Ottawa, Canada 


[Received July 27, 1955] 


THE rate of diffusion of *4Na in NaCl has been measured (Mapother, Crooks 
and Maurer 1950) by evaporating a film of 24NaCl on to the surface of a 
crystal of NaCl and subsequently determining the concentration of 24Na 
as a function of the distance from the surface. Two diffusion processes 
have been found—diffusion by vacancies created thermally at higher 
temperatures (intrinsic range) and by vacancies produced by positive 
divalent impurities at lower temperatures (structure sensitive range). 
Some problems of interpretation raised by the experiments have been 
discussed by the above authors and by Seitz (1954). The study of the 
diffusion of the chloride ion using the same technique with the isotope 
*6C] (Chemla 1952) has been limited to temperatures exceeding 650°c by the 
low mobility of the chloride ion. The present note describes a different 
tracer method which allows diffusion of the chloride ion to be studied at 
much lower temperatures. Preliminary results for the diffusion co- 
efficient and activation energy of chloride ion diffusion at 510°c and below 
are compared with data from the above mentioned papers. 

The method is based on the observation that a rapid isotopic exchange 
takes place between chlorine gas and the surface of small NaCl particles 
(specific surface +30 m?/gm) even at room temperature. A _ similar 
observation on the exchange between HCl and KCl was reported by 
Clusius and Haimerl (1942). It is unlikely that chlorine will dissolve in 
NaCl at temperatures considered here (Mollwo 1937). Therefore, if par- 
ticles of NaCl of any. size containing **Cl are surrounded by non-radio- 
active chlorine gas, the °°Cl diffuses to the surface and appears in the gas 
due to the rapid exchange. Assuming that the concentration of *°Cl at 
the surface is zero (i.e. that the rate of exchange is rapid compared with 
the rate of diffusion) and that the particles are semi-infinite, the number 

_of 36C] atoms in the gas at time ¢ is given by 


N(t)=2AC,»/(Di/7) (Carslaw and Jaeger 1947)... (1) 


where A is the area exposed to the chlorine, C, is the original 
number of 36Cl atoms/cm? in the NaCl and the diffusion coefficient 
D=D, exp(—E/kT). The condition that C=0 at the surface is only 
strictly true when the number of chlorine atoms in the gas is very large 
compared with that in the solid. For small amounts of exchange both 
conditions required by (1) are fulfilled experimentally. For large amounts 
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(or small ratios of gas to solid) an appropriate treatment can be given 
following the method of March and Weaver (1928). 

If a portion of the gas is scanned by a Geiger counter, a counting rate 
R proportional to N is observed. The following ‘ calibration * makes it 
unnecessary to know explicitly the counter efficiency, the actual volume 
of gas scanned by the counter, or Cy. HCl, prepared from a known mass m 
of NaCl containing a concentration Cy of *°Cl, is condensed into the system 
where it produces a counting rate Rj. It may then be seen using (1) that 
the counting rate R is given by 


VIG b on. 6 


where p is the density of the NaCl. Obviously the calibration requires a 
quantitative separation of the HCl. The method used—treatment of 
NaCl with concentrated H,SO,—gave consistent yields of HCl in excess 
of 98%. 

The apparatus used for the exchange experiments requires only a brief 
description. The detailed results given below were obtained for NaCl 
(Merck reagent grade) which had been fused in a platinum crucible and 
to which had been added a small amount of NaCl containing the radio- 
isotope 26Cl. Samples for the diffusion experiments (approximately 0-2 gm 
each) consisting of pieces about 0-2 mm in dimension were sealed in a 
Pyrex apparatus with about 50 cm*s.t.p. of chlorme which had been 


twice distilled from trap to trap. Before sealing, the apparatus was. 


evacuated to a pressure of 1x 10-6 mm Hg. A part of the sealed system 
was the jacketed volume of a Geiger counter and the chlorine was circu- 
lated through this volume by convection. The level of activity was read 
from a counting rate meter connected to a chart recorder. The NaCl 
samples were contained in a Pyrex bucket in a part of the apparatus 
which could be heated by an external furnace. Temperatures were deter- 
mined with a chrome-lalumel thermocouple placed on the wall of the part 
of the apparatus inside the furnace. 

The results are given in figs. 1 and 2. In one experiment (fig. 1) the 
NaCl was held successively at 358°, 435° and 495° and in the second (fig. 2) 
at 467° and 510°. RR? has been plotted against ¢ for both sets of results 
and the linear relation required by eqn. (2) has been found although C 
and D both show a break after 10 to 15 hours. 

The logarithm of the various slopes is plotted against the reciprocal of 
the absolute temperature in fig. 3. It is evident that the initial reaction 
at all temperatures is the same and the activation energy H derived from 
the slope of the solid line is 1-67+-0-05 ev. Within the experimental un- 
certainty the later stages of C’ and D give the same activation energy 
(dashed line of fig. 3). It must be concluded, therefore, that a different 
diffusion process such as diffusion down grain boundaries is not pro- 
ceeding here, but rather a slight change of Dy has taken place. Various 
physical explanations of the break have been considered, viz. lack of 
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The diffusion of *6Cl out of NaCl at 358°c(A), 435°c(B) and 495°c(C). 
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equilibrium of thermally created vacancies, evaporation of NaCl from the 
surface at high temperatures, recrystallization of the surface, but none 


has been fully satisfactory. 
Fig. 3 
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Activation energy plot. A, B, C, D and E as in figs. 1 and 2. 


Using the geometrical area of the NaCl samples (25-+-6 em?) with eqn. 
(2) the solid and dashed lines of fig. 3 lead to the values 0-10 and 0-03 
cm?/sec for Dj. The effective area of the samples as regards diffusion 
might be larger due to grain structure, but could hardly be larger by 
orders of magnitude. The present results are compared with those of 
Chemla (1952) for Cl~ diffusion, and those of Mapother, Crooks and 
Maurer (1950) for Na* diffusion in the table. 


Chloride Ion Sodium Jon 
Present 
work! Chemla?} M—C—M1 M-C-M2 
D,)—cm?/sec 0-10 3x 104 | 1-6x 10-8 3 
H—electron volts 1-67 2-70 0:77 1-80 
H—electron volts 1-03 1-03 
(Dalai las); 3x 105 2x 106 


1Structure sensitive range. 
* Intrinsic range. 
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The low values of H and D, found in the present experiments places the 
diffusion process in the structure sensitive range where the concentration 
of vacancies no longer depends upon the temperature. Hence, the acti- 
vation energy is presumed to be that for diffusion of a negative vacancy 
or of a complex whose rate of diffusion is determined by the diffusion of 
vacancies as assumed by Mapother et al. (1950) for positive ion diffusion. 
The diffusion coefficient may be given more explicitly by 


D=(Do)2 exp (—E/kT)=v_d? exp [—(U+ W/2)k7) (intrinsic range) 
D=(Dg)1 exp (—E/kT)=v_C_d? exp (— U/kT) (structure-sensitive 
range) j 


’ 


(3) 


for the negative ion with similar expressions for the positive ion. Here 
v_is the frequency of vibration of a negative ion, d the interionic distance, 
U the activation energy for motion of a vacancy and W the energy of 
formation of a pair of separated vacancies. C_ is the proportion of lattice 
sites which are vacant in the structure sensitive region. 

AE in the table is, therefore W/2 and, as should be expected, the same 
value results from chloride and sodium ion diffusion experiments. The 
ratios of the diffusion coefficients, (D),/(Do);, are also roughly the same 
for both types of sites. Since (D,)o/(D );=1/C_ for negative ions and 
(Do) 2/(Po)s=1/C+ for positive ions, then the concentrations of negative 
and positive ion vacancies are approximately equal in the structure sensi- 
tive region. According to the expressions (3), (Do)+/(Do)-=v+/v— m the 
intrinsic range. For NaCl, v,/v_-=10-4, which may be compared with 
10-6 for KCl (Seitz 1954) and 10-? for NaBr (Schamp and Katz 1954, 
Mapother, Crooks and. Maurer 1950). The large difference between v, 
and v_ may be expressed in another way using the transition state theory : 

(Do). (AS,— AS_) 


ec ae eae Te ae (4) 


where AS, and JAS_ are the entropies of activation for the motion 
of positive and negative ions. The ratio v,/v_=10™~ leads to 
AS_—AS.,=18-3 cal/mol deg. Possibly this ditference may be accounted 
for by the different polarizabilities of the ions. 

Additional information may be had from the transference numbers of 
Cl- and Na* measured by Tubandt et al. (1931) and quoted by Jacobs and 
Tompkins (1952). These presumably are for the structure-sensitive 
range. The ratio of the transference numbers gives 


t. (AS_—AS,) (Ue U) 
i, =exp [— OxD Sa a ee * 


For Tubandt’s results U_— U_,=0-72 ev as compared with 0-90 ev from 
the data of the table. Furthermore, (4S_—A/S_) is approximately 
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13-7 cal/mol deg. which agrees reasonably well with the value found from 
the data in the table. 

As has been pointed out, the cause of the structure sensitive region for 
positive ion diffusion lies in the presence of an extremely small concen- 
tration of divalent positive impurities (Mapother, Crooks and Maurer 
1950). The foregoing discussion suggests that these might also affect the 
negative ion diffusion. However, purification of the NaCl by growth of 
erystals from the melt led to no significant change in # and Ds for the 
structure sensitive range. Spectrographic analysis of this recrystallized 
NaCl showed a content of less than 0-005°%, of each of nine metallic im- 
purities. Further, the addition of 0-3°, of Cd++ ions (in the form of CdCl,) 
produced no appreciable change in H and Dy. One is obliged to say, 
therefore, that chloride ion diffusion in the structure sensitive region 
proceeds by a mechanism independent of positive ion impurities. Eluci- 
dation of the mechanism must await further experiments. The possible 
existence of multivalent negative ions (e.g. O-~, S—_) in the lattice has not 
been investigated yet. In analogy with positive ion diffusion, it might 
be anticipated that multivalent negative ion impurities could affect the 
chloride ion diffusion. 

Finally, brief mention should be made of the exchange between chlorine 
gas and the small particles of NaCl. The particles were prepared as des- 
cribed by Thomson ef al. (1955) and from electron micrographs appeared 
to be cubes presumably bounded by {100} faces. At room temperature 
a rapid exchange reaction involving about 50°, of the surface chloride ions 
was observed. A much slower, apparently diffusive reaction followed. 
Between room temperature and 80°c the activation energy for the slow 
process was about 0-3ev. The measurements at higher temperatures 
were complicated by the occurrence of sintering and hence estimates of D, 
were quite unreliable. ‘The low activation energy suggests surface rather 
than bulk diffusion. The fact that not all of the surface chloride ions 
entered into the rapid exchange is to be expected since perfect {100} faces 
would exchange only very slowly with adsorbed chlorine molecules. The 
rapidly exchanging portion of the surface may then be associated with 
disorder of the surface producing ‘exposed ’ NaCl molecules which are 
easily exchangeable. 
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ABSTRACT 


Experiments on single crystals of cadmium at temperatures between 
1-2° and 4-2°K have shown that measurable transient creep occurs, and 
that the amount of such creep does not fall off rapidly with decreasing 
temperature. It is suggested that this creep is due to the quantum 
mechanical tunnel effect, which allows dislocations to penetrate barriers 
to their motion that cannot be overcome by the action of the applied 
stress alone. 


§ 1. [XYTRODUCTION 

‘Iv is now 25 years since Meissner, Polanyi and Schmid (1930) first 
reported results of tests on metal single crystals immersed in liquid 
helium. Their main concern at that time was to show that plastic 
deformation still occurred at stresses well below the ‘ theoretical’ yield 
stress for a perfect crystal at temperatures very close to absolute zero, 
and that plasticity was therefore not a purely thermal effect. Since 
that time dislocation theory has provided an adequate explanation for 
the discrepancy between the observed and the * theoretical ’ yield stress, 
but one of their results is of interest in a rather different connection. 
They noticed that, even at these temperatures, there seemed to be a. 
small amount of transient creep; that is to say that, after application 
of the stress, the strain continued to increase for some time. In their 
paper they showed two creep curves for cadmium crystals, one at 4:2°K 
the other at 1-2°K and the creep in these two cases was not very different 
in magnitude. 

Recent theories of transient creep all assume that the phenomenon is 
due to thermal fluctuations assisting the applied stress to force dislocations 
past some obstacle, and therefore on these theories creep should vanish 
as the temperature approaches absolute zero, unless zero point energy 
can supply the necessary fluctuations in the same way as thermal energy. 
As this point is of some interest in the theory of creep, and as no further 
experiments seem to have been reported to confirm the results of Meissner, 
Polanyi and Schmid (which they themselves admit might have been due 
to accidental fluctuations of the stressing system), an attempt has been 
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made to investigate this phenomenon, and to find out somewhat more 


about this low temperature creep. This paper presents a preliminary 
account of the results. : 


Fig. 1 


Schematic diagram of the testing machine used. 


§ 2. SPECIMENS AND APPARATUS 
The crystals used were made of Johnson Matthey spec-pure cadmium 
(batch no. 4753X) by the Andrade technique. Lengths of wire about 
10 cm long and 2 mm diameter were sealed in an argon atmosphere in 
glass tubes, of which the diameter was somewhat greater than the wire 
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diameter, and a travelling furnace was then passed over the tubes to 
melt the wire and allow it to resolidify from one end. The resulting 
crystals were long enough to allow two or more short specimens to be 
cut from each. The cutting was done with a flame, and the specimens 
were melted onto copper loops at each end. Back reflection Laue 
patterns were used to determine the orientation. 

The apparatus used was a simple testing machine, specially designed 
for operation in liquid helium. A diagrammatic sketch of the assembly 
is shown in fig. 1. The load from the spring A, which could be extended 
using the winch B, was applied to the specimen S through the thin 
German silver tube C. The movement of the top of the tension system 
was measured directly by means of the telescope T, which was focused 
through a window onto an electron microscope graticule D fastened 
onto the wire which transferred the load from the spring A to the tube C. 
The telescope could be moved up and down with a micrometer, whose 
reading gave the movement of D to an accuracy of 0-002 mm. A second 
graticule, shown schematically at E, was also measured to allow for 
bodily movements of the machine relative to the telescope. The bottom 
of the specimen was hooked to the small threaded rod F, which carried 
a nut that could be screwed up to the bar G at the beginning of a test. 
G was connected by two rods HH to a second German silver tube [, 
which was connected at its top end to the main case J. This case was 
airtight, so that, when the specimen had been mounted and the vacuum 
flask K put in position with a rubber band L, the whole assembly could 
contain helium. A second vacuum flask containing liquid oxygen was 
placed round the helium flask. A stop M was provided on the tension 
wire, so that, if the specimen broke, it would come against the adjustable 
annulus N. The whole assembly, other than the telescope, was removed 
bodily to the helium liquifier for filling with liquid helium. 

The design of the machine was intentionally kept as simple as possible, 
so that modifications could be made after the magnitude of the creep 
was known, and for this reason there are some disadvantages. The 
extension of the specimen shortens the spring, and therefore reduces 
the load, and this not only involves a calculation to deduce the load, 
but also means that a slight fall in load occurs during creep, although 
with the small amounts of creep expected this effect is not large. Also 
the long distance between the measuring point and the specimen means 
that the elastic deformations of the machine, determined from a blank 
test, must be taken into account when deducing the stress-strain curve 
from the data. However, as neither of these points stop the machine 
from demonstrating the existence of creep, it was decided to make a series 
of experiments without further modification. 


§ 3. RESULTS 


Four specimens have so far been tested, some on more than one occasion. 
One of these was made of commercial purity cadmium, and the other 
three from the spec-pure material. In all cases there was no doubt that 


Creep of Cadmium Crystals at Liquid Helium Temperatures 403 


creep occurred. When the winch was moved to its next stop, the graticule 
moved up a large distance, and also oscillated sideways with the rod. 
but if the cross-wires were set on it, it could be seen to continue rising, 
and after the oscillations had died down, it continued slowly to rise for 
the next half-hour. As it seems unlikely that stray oscillations picked up 
from the laboratory floor could have been of the magnitude of the initial 
oscillations, the continued stretching after the initial oscillations were 
over is strong evidence for the existence of true creep, not caused by the 
load attaining a new maximum in oscillations. 

Figure 2 shows the load-extension curves obtained. The amount of 
creep observed at any one setting of the winch is marked by the roughly 
horizontal lines. Specimens 1 and 2 were cut from the same crystal 
(x=38°, A=39°), specimen 3 from another, whose Laue pattern showed 
pronounced asterism, while specimen I was made from commercial 
purity cadmium. 

A number of the creep -curves are shown in fig. 3. It will be noted 
that the creep falls off very rapidly, and in fact appropriate plots showed 
that, while the strain is more nearly proportional to In¢ than to #1, 
it actually fell off faster than predicted by either of these laws. This 
could be attributed to the drop in load as the specimen extends, as this 
drop, though small (of order $°), is greater than that needed to keep 
the stress constant, but similar creep curves at the temperature of liquid 
nitrogen seemed to be more nearly proportional to In ¢, and reasons 
will be advanced below for expecting such a rapid decrease theoretically. 
It is also interesting that decreasing the temperature during a test did 
not affect the amount of creep by anything like the ratio of the 
temperatures. Unfortunately, pumping off helium vapour to reduce the 
temperature caused severe oscillations, both from the pump itself, and 
also from the violent agitation of the helium until the A point was reached. 
For this reason the stress was usually lowered somewhat before pumping 
was started, and only raised again when the specimen was down to the 
lower temperature and the helium was still. In all cases the load—extension 
curve appeared to be raised above its extrapolated course before cooling. 
so that some temperature-dependent effect remains even at these 
temperatures. A further difficulty arose in these tests below the A point : 
the whole machine tended to move due to forces in the connection to the 
pump, and so the changes in the zero were more serious than in the tests 
at 4:2°x. This further reduced the accuracy in these creep curves, a 
fault that it is hoped to remedy in further experiments. 

A few tests were also performed at liquid nitrogen or liquid oxygen 
temperature to find the magnitude of the creep at these temperatures. 
The spring used was rather harder than was really desirable for tests on 
single crystals, but the creep found was very nearly proportional to 
Int, which further shows that the failure of this law at helium temperatures 
is not due to the machine. As yet tests at intermediate temperatures 


have not been possible. 
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Load-extension curves: (a) shows the curves for tests 1 and 2, which were 
on specimens cut from the same crystal, of which the initial orientation 
was x¥=38°, A=39°. (b) shows the results of test 3 in which a specimen 
cut from a crystal which showed pronounced asterism, was tested 
on two different occasions with a room temperature anneal between. 
(c) shows the curves for specimen I, made from commercial purity 


cadmium, and tested at various temperatures with room temperature 
anneals between each. 
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Some of the creep curves obtained. The curves from each test are shown on 
a separate diagram. (d) Test 1, (b) test 2, (c) test 3, (d) test 3a, the 
same specimen as test 3 re-tested after a room temperature anneal, 
(e) test I, (f) test I, creep curve obtained at room temperature. Note 
the different scales in this case. 
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§ 4. DiscussIon 


The first question that must be answered is whether the creep observed 
might not be the normal, thermally activated creep. Wyatt (1953) has 
found experimentally that at low enough temperatures (usually about 
the liquid nitrogen range) most metals obey a creep law of the form 


«=a In (yt+1), ny ea 


where « is proportional to the absolute temperature, and it is difficult 
to see how, on any theory depending on thermal activation for the creep, 
the temperature variation could be slower than this. Furthermore, a 
theoretical estimate of the value of « can be obtained from the theory 
of Mott (1953). On this theory 

kT /\ do 
On = Pp / a dé 3 ° . . 5 . * (2) 


where B is the activation energy needed to overcome the elementary 
barrier to flow, and this can be expected to be about ley, both on 
general grounds, and also in agreement with creep measurements at 
higher temperatures ; do/de is the gradient of the stress-strain curve. 

In order to compare the present results with these predictions, it is 
convenient to deduce a value of « from the difference between the strains 
at two times ¢, and ¢t, after the application of the load. Using eqn. (1) 
above 


yto+1 
9— €4=« In— , 3 
€o—€4=4 Del (3) 
and if yt,>1 
Come Gia In (ty/t,). . . . . ° . (4) 


Furthermore, if yt,;> 1 ‘is not true’, the value of « obtained from ,(4) 
will be smaller than the correct value obtained from (3). 

The values of « obtained using (4) for various creep tests are shown 
in the table, together with the values expected on Mott’s theory if B is 
taken as lev. It will be noticed that only for temperatures above the 
helium range is the value of « of the same order of magnitude as the 
expected values apy. At the low temperatures the value of « is too 
large by about a factor 10, and is more nearly independent of temperature 
than proportional to it. 

As the usual theories of creep fail to account for the creep at helium 
temperatures, it is necessary to consider what effects have been neglected 
in their formulation, and one of these is zero point energy. It is not, 
however, sufficient simply to add some term, such as k@ to the term k7, 
as creep presumably depends on the cooperative movement of many 
atoms, and the spectrum of zero point energy is very different from that 
of the thermal oscillations. It is perhaps more enlightening to consider 
the effect, not as due to the zero point oscillations of the lattice, but as 
due to the tunnel effect through the energy barrier trapping the 


Ea 


——— 
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Results of Creep Tests Analysed to show the Amount of Agreement 
with the Theory of Temperature Activated Creep 


Test Temp. Load peo ty te o oTH | - 
(°K) (g) ade | (sec) | (sec) | (105) | (x 105) 
Spec-pure moe 
1 4-2 480 6-6 1 60 44 : 
1 4-2 538 6-2 WN 1s 53 a 
l 4-2 715 eee dat 60 43 6-8 
1 77-4 393 18-2 1 60 37 36-7 
Ye 4-2 425 6-6 il 30 Pe 5-4 
7 4-2 480 5-4 1 30 40, 6-6 
2. 4-2 544 7-4 1 42 60 4-8 
2 1-4 714 4-7 1 43 24 2-5 
2 1-4 798 rei a 15 30 127 
Oy. 1-4 861 4-5 0-5 7:5 il il 
33 4-2 620 3°7 il 82 40 9-8 
3 4-2 677 3°5 i 19 330) 10-4 
3 4-2 734 O25 1 PA 44 11-0 
3 1-2 920 2°8 0-5 14 30 4-2 
3 1-2 983 ay 1 46 65 4-4. 
3a 4.2 442 4-] if 1s 12 8:8 
3a 4-2 490 |- 3-8 1 24 26 9-5 
3a 1-2 617 4-3 1 34 20 2-4. 
Commercial 
purity 
I 4-2 1268 18-5 1 50 29 1-9 
I 4-2 1437 15-6 6 74 36 2°3 
I 4-2 1519 14-4 1 60 28 2-5 
I 90 1014 13-1 0:5 56 40 59 
I 90 1176 11-5 il 93 66 67 
Steady 
I 300 620 == = || Oo 6 783 state 
creep 
I 90 845 30 1 3 49 25-8 
I 90 927 18-2 0-5 4 41 43 
I 90 1013 13-7 0-5 4 yap} 57 
I 90 1092 12-1 0-5 9 36 64 
I 90 1170 11-4 0-5 3 70 68 
i 90 Wot 10-7 0-5 4 70 72 
I 90 1330 10-1 0-5 4 62 76 
ie F920 1490 7:8 0-5 60 64 100 
I 90 1570 7-4 1 pag 61 104 


dislocation. The dislocations that are responsible for the first creey 
are those for which the applied stress is very nearly adequate to take 
them over their barriers, and therefore the barriers facing them are very 
low and very narrow. To make a proper calculation of the creep which 
can arise from this effect would involve, among other things, an exact 
knowledge of the shape of the potential barrier faced by dislocations 


408 Creep of Cadmium Crystals at Liquid Helium Temperatures 


and of the effective mass of dislocations, but it is certain that some 
dislocations, unable to surmount their potential barriers with which they 
are faced with the aid of stress alone, will be able to penetrate them, 
given time, by the normal quantum mechanical tunnel effect. The 
oscillations of the lattice in modes other than that considered as the 
movement of the dislocation will probably assist the tunnel effect, as 
the barrier will itself fluctuate in height due to their action. The only 
question is whether this effect produces enough creep to be observable. 
While the theoretical solution of this problem must be left for the 
present, the experiments described above seem to indicate that some 
effect such as the tunnel effect can give measurable creep. The creep 
produced in this way should be independent of temperature, and it 
should also fall off with time more rapidly than would be the case were 
the same dislocations overcoming the same barriers with the aid of 
thermal motions, since a higher barrier is also broader, and the tunnel 
effect falls off exponentially with both height and width. In both these 
respects the observed creep agrees qualitatively with the predictions. 
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ABSTRACT 


X-ray measurements are reported of the thermal expansion of 
industrial and gem-quality diamond from —150°c to +950°c. It is 
found that the thermal expansion of industrial diamond shows an anomaly 
at low temperatures. No such anomaly is found for gem-quality diamond, 
for which Griineisen’s law holds in general throughout the range. 
The values of the unit-cell dimension at 20°c are: industrial diamond, 
y=3-5669, A, gem-quality diamond, a,=3-56684 A. 


§ 1. INTRODUCTION 


THE value of the thermal expansion coefficient of diamond and its variation 
with temperature are of considerable interest in the study of the graphite— 
diamond equilibrium diagram since the thermal expansion is directly 
related to the difference between the specific heats at constant pressure 
and constant volume. Berman and Simon (1955) and Liljeblad (1955) 
have recently put forward suggested diagrams ; and Berman and Thewlis 
(1955) have shown that Berman and Simon’s diagram is to be preferred 
in the light of recent measurements of thermal expansion up to 900°C. 
It is the purpose of the present note to describe these measurements. 


§ 2. MrerHop 


X-ray powder photographs were used in making the measurements, 
the temperature range covered being from —150°c to +950°Cc, a greater 
range than has previously been achieved in one series of experiments. 
The low-temperature photographs were taken with a modified 19 cm 
Unicam camera (Thewlis and Davey 1955) in which the specimen is 
cooled by a stream of cold nitrogen, and the high-temperature photo- 
graphs with a Unicam high-temperature camera of early design, of which 
the specimen holder had been modified. . 

Two batches of diamond were examined, industrial and gem-quality, 
the latter being ground to give crystallites less than 5 microns in diameter. 
The material to be examined was held in thin-walled silica capillaries of 
internal diameter 0-5 mm. 
ec NR ete oe tan sR kd SE aI ae OY ee 
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The approximate temperature was measured in the low-temperature 
camera by chromel—alumel thermocouples and in the high-temperature 
camera by platinum—platinum—rhodium thermocouples. The actual 
temperature was measured, however, by measuring the lattice parameter 
of silver, pure silver powder being mixed with the diamond powder so that 
both sets of reflections occurred on the same photograph. For tempera- 
tures below room temperature the lattice parameter values used were 
those given for silver by Owen and Williams (1954) ; above room tempera- 
ture those of Hume-Rothery and Reynolds (1938) were used. Un- 
filtered Cu K radiation was employed throughout except for a few high- 
temperature photographs of industrial diamond, where CoK was 
substituted. 

The accurate value of the unit-cell side for each photograph was deter- 
mined by plotting the observed value against the Nelson—Riley function 
(Nelson and Riley 1945) for a number of reflections and extrapolating to 
6=90°. Independent extrapolations for the Ka,, Ka, and Kf reflections 
were carried out to allow for systematic errors in reading the positions of 
the K«, and Kx, reflections (Thewlis 1955), the mean value being taken 


as correct. 
§ 3. RESULTS 


Figures 1 and 2 show the values of the unit-cell dimensions at different 
temperatures for industrial and gem-quality diamond repectively. It 
wlll be seen that the high-temperature values are very similar but that, 
for industrial diamond, there is an anomaly at about —30°c. This 
anomaly is brought out more clearly in fig. 3, where the linear coefficient 
of expansion, « (i.e. (1/a)(da/dt) where a is the unit cell side and ¢ the 
temperature), is plotted against temperature for industrial diamond. 
None is present in fig. 4, which shows a plot of « against temperature for 
gem-quality diamond. 

It will be seen, as is to be expected from figs. 1 and 2, that the two types 
of diamond give very similar values of « at the high temperatures which 
are of interest in connection with the graphite-diamond equilibrium 
diagram, but there is a small region near 0°c where, for industrial 
diamond, the coefficient of expansion becomes negativet. The anomaly 
could presumably arise from some intra-crystalline impurity, and it is 
interesting to note that previous hints of anomalies in this temperature 
region are not unknown. For example Réntgen (1912) refers to a density 
maximum at —41-7°c, the observation of which is attributed to Fizeau. 


§ 4. Discussion 


Figure 2 shows, as well as the results of the present measurements, the 
values obtained by Krishnan (1946) and Straumanis (1951), the only 
other workers, as far as we know, who have measured the unit-cell 
dimensions of diamond at different temperatures. It will be seen that 
the agreement is satisfactory. 


+The values of « in the anomalous range are not precise, and the slight 
negative value should not be regarded as established beyond doubt. ; 
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In fig. 4 the results of these and other observations of the thermal 
expansion coefficient of diamond are shown. Joly’s (1897) results differ 
widely above 400°c from ours but consideration of his method (the 
measurement of the optical projection of a small heated diamond) shows 
that this is hardly surprising. Krishnan’s values, and those of 
Dembowska (quoted by Réntgen 1912), are in fair agreement with ours, 
although the shoulder at about 250°c is less pronounced in Krishnan’s 


Fig. 3 
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curve. The results of Fizeau (1869) and Straumanis correspond to such 
restricted temperature ranges and to such relatively few actual measure- 
ments that little comment is called for. It is our opinion that the true 
values for the linear expansion coefficient of diamond up to 600°C lie 
between the Krishnan—Dembowska curve and ours. Above that 
temperature the consistency of the observed values of the unit-cell 
dimensions, as shown in fig. 2, suggests that the errors in our curve are 
small. Berman and Thewlis (1955) point out that this curve agrees with 
what would be expected if Griineisen’s law held, and is the justification 
for extrapolating to higher temperatures according to that law, 


& (x10°) —> 
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§ 5. Nore on THE ACCURACY OF MEASUREMENT OF THE 
Untt-Ceti SIpE 
As will be seen from figs. 1 and 2 the random errors of individual 
measurements of the unit-cell side do not exceed 0-0001 4, and are 
usually less. This is consistent with the degree of accuracy attained in 
our previous work (e.g. Thewlis 1955). Our value for the unit-cell side of 
gem-quality diamond at 20°c is 3-5668, 4, the literature values ranging 
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from 3-5667, A (Straumanis) to 3-5669, A (Lonsdale 1944). For the unit- 
cell side of industrial diamond at 20°c we obtain the value 3-5669, A. 

No allowance has been made for refraction, as the correction on this 
account is only 0-00001 A. The wave-length of Cu Kz, radiation has been 
taken as 1-54050 A. 
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Clarendon Laboratory, Oxford} 


[Received December 24, 1955] 


Borvont (1949 and 1954) has reported the existence at low temperatures 
of a large peak in the values of internal friction plotted against temperature 
for several cold worked metals, and we are making a more detailed survey 
of the position. Our work has so far been confined to experiments on 
99-999°% “polycrystalline copper, and it has been shown (Niblett and 
Wilks 1955) that the different types of results obtained by Bordoni for 
different metals can all be induced in copper by subjecting it to varying 
amounts of cold work. Although measurements are still in progress it 
seems worth while to give a preliminary report, particularly of those results 
which appear to be inconsistent with an explanation of the peaks, published 
recently by Mason (1955). 

All the measurements have been made on bars of copper which were 
first annealed in argon for one hour at 600°c and then plastically deformed 
in a tensometer. The ends of the bar, which had been gripped in the jaws 
of the tensometer, were cut off and the bar mounted in the cryostat on the 
edges of razor blades. Measurements of the internal friction were then 
made by driving the bar electrostatically into transverse vibration at its 
resonant frequency (about 1100 c/s) and observing the subsequent decay of 
the oscillations. The strain amplitudes, that is the strain due to the 
vibrations, have always been small ; of the order of 10~*. All our results 
have been corrected for damping due to the thermo-elastic effect, but this 
is in any case small. 

It turns out that the friction depends on the prior cold work in a rather 
complicated manner, but for deformations greater than about 2°% a curve 
of the type shown in fig. 1 is obtained. This behaviour is similar to that 
reported by Bordoni for copper except that his measurements do not seem 
to have been sufficiently sensitive to show the subsidiary peak. For lesser 
amounts of cold work the two peaks are still observed, but are smaller 
and appear to be superposed on a background of friction arising from some 
obviously different mechanism (fig. 2). This background friction is also 
of interest for its own sake and further investigations are in progress 
particularly to determine more fully the complex amplitude dependence. 

The peaks appear to be of the type associated with relaxation 
phenomena. Thus the friction is to a first approximation independent of 
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strain amplitude, while the positions of the maxima depend on the 
frequency. A thinner bar vibrating at about 400¢/s had its main 
maximum some 7° lower than that shown in fig. 1 ; such ashift corresponds 
to an activation energy of 0-07 ey and this value is also consistent with the 
shape of the friction versus temperature curve. (However, if one should 
also use Bordoni’s peak in copper at 30 ke/s to determine the energy one 
arrives at a somewhat higher value, say about 0-1 ev.) As the specimens 
were deformed at room temperature and left there for several days prior 
to an experiment, this activation energy cannot correspond to the motion 
of a free defect such as a vacancy or interstitial for in this case the defects 
would have diffused out. Thus the most likely source of the damping is 
the motion of a short segment of dislocation line, as is confirmed by the 
fact that the friction does not anneal out save by heating to 350°c (see 
fig.-1). 
Fig. 3 


fe) 50 100 150 200 


Temperature °K 


The internal friction of copper, containing 0-0026% Bismuth, and 0-032% 
Phosphorus, strained 5-5°%. 


A detailed mechanism has been proposed by Mason (1955) which, he 
claims, gives rise to an activation energy of the right order of magnitude. 
He envisages the dislocations as being anchored at pinning points, and that 
the lengths between these points suffer small displacements. Hence he 
arrives at an activation energy which is proportional to the distance / 
between the pinning points for the length of dislocation is question. As 
this length / can presumably take up a very wide range of values, it is not 
clear how such a mechanism can give rise to a fairly well defined activation 
energy, with a characteristic peak at a given temperature. Moreover, if the 
pinning is due to the interaction of dislocations then the distance / should 
decrease with increasing cold work, but as mentioned above the shape and 
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position of the peak in our experiments is virtually independent of the 
amount of prior deformation. Alternatively, if the pinning is due to 
impurities then the addition of further impurities should reduce the value 
of / and give rise to a lower activation energy. Our results for copper 
containing appreciable amounts of phosphorus and bismuth and strained 
5-5% are shown in fig. 3. Although the amplitude of the peak is consider- 
ably reduced, the shape and position of the curve is substantially the same 
as that for pure copper, if allowance is made for the relatively larger 
contribution of the background damping. Thus the present results are not 
in agreement with Mason’s treatment. 

A much more promising explanation has been given by Seeger (1955) who 
considers the displacement of a dislocation line lying parallel to the lattice 
lines. He has shown that the initial motion will consist of the displace- 
ment of a short segment of line whose length is such that the resulting 
energy change isa minimum. Thus the vibrations of the bar set in motion 
sections of dislocations of a unique length and with a unique activation 
energy, and this energy may be shown to be of the right order of magnitude. 
Moreover, as in a f.c.c. lattice there are two principal directions in which 
the displacements can take place, we would expect there to be two peaks 
rather than one. 

Finally we note that while the increase of the friction with increasing 
cold work is almost certainly due to the increasing density of dislocations, 
some other mechanism is needed to account for the fact that the maximum 
friction observed after about 8°, deformation is rather less than that after 
2%. Of course, a decrease of internal friction is often observed following 
largish amounts of cold work and is usually, if somewhat vaguely, 
accounted for by supposing that the dislocations become tangled up with 
one another. In the present case however it is clear that the interactions 
which reduce the friction must be such as not to affect the value of the 
activation energy. Further experiments to elucidate the nature of this 
mechanism are now in progress. 
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ABSTRACT 


Measurements have been made of the magnetic polarization of 450 
specimens of basaltic lavas of the Deccan Trap taken from two sites 
about 500 miles apart. With the exception of some nearly randomly 
magnetized specimens, the main direction of magnetization is N. 155° E. 
and the dip is 53° downwards. The most plausible interpretation of 
this result is (a) that India has drifted north from a position about 34° 
south of the equator when the rocks were formed some 70 million years 
- ago and has rotated anti-clockwise through 25°, and (b) that either the 
earth’s field was reversed when the rocks were formed or that the rocks 
became magnetized in the opposite direction to the field by some physico- 
chemical mechanism. 

The pole position corresponding to these results does not agree with that 
found by other workers for British and American rocks of the same period. 
This suggests that a movement of India, relative to North America and 
Europe, has taken place at some time during the past 70 million years. 


§1. LyTRODUCTION 


In October 1954 the Physics Department of Imperial College embarked 
on a three-year project of palaecomagnetic research in India. The 
scheme arose as a result of a previous survey made in Britain when it 
was found that certain sediments laid down over a wide area during the 
late Palaeozoic and early Mesozoic eras are consistently magnetized along a 
north-east south-west axis, and had a considerably lower inclination 
than the present earth’s field (Clegg, Almond and Stubbs 1954). This 
result could only be explained with any degree of plausibility on the 
assumption that the land mass of England has drifted northward, and 
has rotated approximately 30° in a clockwise sense, relative to the 
earth’s geographical axis, at some time since the end of the Triassic 
period. Further evidence of past movements of the British land mass has 
also been obtained from palaeomagnetic measurements made by Creer, 


Irving and Runcorn (1954), and the data now to hand are sufficiently 
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conclusive to encourage the belief that by carrying out palaeomagnetic 
surveys on a global basis it may-be possible to determine the positions 
occupied by the various land areas of the earth relative to the poles during 
different geological epochs. Recent work on American rocks by Du Bois 
(1955), Runcorn (1955) and Graham (1955) lends strong support to this 
expectation. 

Among the areas of the world especially suitable for rock magnetic 
studies is the sub-continent of India. Of especial significance is the 
permo-carboniferous glaciation which suggests that the geographical 
latitude of India may not always have been what it is to-day.f Moreover, 
Indian geology has been studied in detail, and numerous exposures exist 
revealing a large part of the stratigraphical column. In particular, the 
Indian peninsular contains the important volcanic formation known as 
the Deccan Trap. This extensive series of basaltic lava flows covers an 
area of some 200 000 square miles, and was most likely extruded during 
the late Cretaceous and early Eocene, and so prior to the culmination of 
the Alpine orogeny when the great tectonic movements responsible for 
the raising of the Himalayas were in progress. Its advantages as a field 
for palaeomagnetic work were first pointed out by Irving (1954), who 
obtained a few samples of the Deccan Trap lavas from different localities, 
and found them to have a consistent axis of magnetization, directed 
roughly south in azimuth, but dipping towards the south. Despite 
the small number of specimens examined, the consistency of the results 
suggested that a more extensive magnetic survey of these rocks would be 
of interest, and we decided, in planning the present expedition, to con- 
centrate initially on the Deccan Trap area. _ 

The first phase of the work has now been completed. Two of us 
(EK. R. D. and D. H. G.) spent two months in India during the winter of 
1954-1955 and collected, in all, 208 rock samples. These were taken from 
a number of lava flows occurring in two different localities some 550 
miles apart ; the first in the neighbourhood of Linga near the District 
Capital of Chhindwara in the Central Provinces, and the second close to 
Khandala on the Bombay—Poona highway. The two areas will be 
considered separately. 

The rock samples, which weighed in all about half a ton, were shipped 
to England for measurement. The standard method used for cutting 
cylindrical samples from the original blocks, and of determining the 
direction of magnetization using a simplified version of the astatic 
magnetometer described by Blackett (1952), has been described elsewhere 
(Clegg, Almond and Stubbs 1954). 


2. RESULTS FROM THE LinGA AREA 
Over an area of some 50 square miles around Linga there are a number 
of isolated exposures of basaltic lavas of the lower Deccan Trap series. 
The district has been mapped by Fermor and Fox (1916), who recognized 


j For summary of some of the palaeo-climatological evidence see Holmes 
(1952) and Gutenberg (1951). 
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five separate flows which they designated by the numbers, 1, 2, 2a, 
3 and 4 in ascending order. They estimated the age of these rocks to 
be late Cretaceous or possibly early Eocene and considered that the 
lowest flow (no. 1) might well be the oldest in the Deccan Trap series. 
A total of 69 samples was taken from nine sites including stream beds, 
quarries and road cuttings, spaced a few miles apart. At the first six 
sites to be visited it was possible to recognize each of the individual 
exposures as belonging to one of the flows 1, 2 or 2a. At the remaining 
three sites flows 3 and 4 were exposed, but it was not easy to differentiate 
between the two. 

Laboratory measurements showed that all the samples with the 
exception of those from flow 2 have preferential directions of magnetization 
towards the south-east with rather steep downward dips. The combined 
results of the measurements of declination and dip for flows 1, 2a, 3 and 
4 are shown as histograms in figs. 1 (a) and 1 (b), and on a Schmidt pro- 
jection in fig. 3 (a). 

The samples from flow 2 were found to have random directions of 
magnetization, and the results have therefore been omitted from the 
histogram, but are included as open circles in the projection. 


§ 3. RESULTS FROM THE KHANDALA AREA 


Close to Khandala, on the steep face of the Bhor Ghat, a large vertical 
sequence of basaltic lava flows are exposed in a relatively small area. 
They are thought to belong to the upper part of the middle section of the 
Deccan Trap, and according to Fermor and Fox (1916) they were probably 
laid down during the early —Hocene period, and are younger than the 
Linga rocks. 139 rock samples were collected from 20 separate sites on 
the hillside, at heights ranging from 475 to 2525 feet above sea level. 
Since the individual exposures were on the whole small and scattered, and 
the flows numerous and often no more than a few feet in thickness, it 
was impossible to identify the latter individually. 

The complete results for the whole Khandala area are shown in figs. 2 (a), 
2(b) and 3(b). It can be seen from the histograms and the projection that 
the majority of the specimens are again magnetized in a general south— 
easterly direction with rather steep downward dips. A significant 
feature of the distribution in fig. 3 (b) is the group of points in the north— 
east quadrant and which produce the small secondary, maximum in 
fig. 2 (a). They were all cut from four blocks, two of which came from 
the lowest flows at the base of the Ghat, and the other two of which were 
collected from the top of the mountain at heights of 2275 and 2525 feet. 
Tt will be seen from fig. 3 (b) that these anomalous specimens are magne- 
tized in a direction considerably closer to the direction of the present earth’s 
field than are the great majority of the specimens. It is therefore possible 
that their polarizations have been pulled round towards the present 
field, owing to partial magnetic instability. There are also three specimens 
from this area which show upward dips. These again were taken from a 
single clock belonging to one of the uppermost flows, 
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§ 4. Discussion oF REsuLts 


Further details of the results for the two areas are given in the table. 
Although the declinations of the individual samples were measured with 
respect to magnetic north, the mean values shown in column 7 of the table 


Fig. 3 
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are expressed in degrees east of true north. In the particular region 
concerned the declination is very small, being 0-6° west. The lava flows 
_ sampled were all very nearly horizontal, and the magnetic dips (column 8) 

were therefore measured relative to a horizontal plane. This marked 
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horizontality over large areas is a characteristic feature of the Deccan 
Trap lavas, which are believed to have been laid down as level flows and 
to have remained free, for the most part, from subsequent folding or 
tilting (Fermor and Fox 1916). 


Fig. 3 
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Schmidt Projection 


Circles and dots, lower hemisphere: @ O; Crosses, upper hemisphere: + ; 
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The figures given in the table show that the general direction of mag- 
netization is similar in both areas. However, the rocks from Khandala 
have, on the whole, a significantly bigger downward dip and a more 
easterly declination than those from Linga. This difference is partly 
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due to the anomalous and probable partially unstable specimens from 
the Khandala, which have already been discussed. If these 23 anomalous 
specimens are omitted, the mean values for declination and dip for 
Khandala obtained from the remaining 210 specimens become 147° 
and +-58° (shown in brackets in the table), which are somewhat closer to 
Linga figures. 

The origin of the random scatter of the specimens from flow 2 in the 
Linga area is unknown. Lightning, weathering or some other disturbance 
are possible causes. 

Though the difference between the results from the two sites is still 
significant and is under further study, for the purposes of the present 
discussion we will use the mean value for the two locations. Thus we 
get for the direction of the magnetic intensity (north seeking pole) 
when the rocks were magnetized, the following figures : 


Declination IN. 155° _E: 
Inclination + 53° down. 


It will be seen from figs. 3 (a) and (b) that the direction of magnetization 
of most of the rocks is almost perpendicular to that of the present magnetic 
field in India. The majority of the rocks, with the possible exceptions 
which have already been mentioned, must therefore have had a 
considerable degree of magnetic stability over long periods of time, so 
that it is likely that their polarizations give a true indication of the 
earth’s field at the time of cooling. The relatively narrow peaks in both 
declination and dip shown in figs. 1(b) and 2(b) strongly suggest 
a high degree of stability for most of the specimens. For when magnetic 
instability in rocks is found it generally affects different specimens 
differently, leading to a wide distribution in measured directions. 

We conclude that throughout the whole period during which the lava 
flows of the Linga and Khandala were being laid down, the earth’s 
magnetic field was directed either along a south-east axis dipping down- 
wards, or, if the rocks were self-reversed, along a north-west direction 
dipping upwards. Since the period of time over which the Deccan Trap 
was extruded was probably much longer than the period of the secular 
variation, we can reasonably assume that the mean direction of magnetiza- 
tion of all the measured flows should be close to that corresponding to an 
axial dipole, which we assume to be co-axial with the geographical axis of 
the earth. 

Our results therefore suggest that the Indian peninsular occupied a 
very different position relative to the geographical poles during the late 
Cretaceous and early Eocene periods from that at the present time. In 
particular, in order to account for the rather steep downward dip of 55° 
compared} with the present 25°, it is necessary to suppose that India 


Ae Sa a es ee ee ee ees 
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is 36°. 


2F2 


428 J. A. Clegg, E. R. Deutsch and D. H. Griffiths on 


was then situated considerably farther from the equator. For if the 
earth’s field corresponded on the average to that of an axial dipole, the 
relationship between the mean magnetic dip, 6, and the geographical 
latitude, A, would be given by 


tan A= tan @. 


Substituting the mean values of @ of 53° from column 8 of the table, we 
obtain the values of 34° for A. After taking into account the direction and 
sign of the dips, we therefore have two possibilities, namely : 

(a) This particular area of Central India may then have been about 
34° south of the equator. In this case it is necessary to assume either 
that the direction of the field was opposite to that at the present time, 
or alternatively that the rocks took on a reverse magnetization on cooling. 
We must therefore suppose that the whole sub-continent has moved 
northward across the equator at some time during the intervening 70 
million years between that time and the present day. The south-easterly 
declinations indicate that this transitional movement was accompanied 
by an anti-clockwise rotation of 25°. 

(b) The observed directions of polarization are also compatible with 
the assumption that the rocks were formed in the northern hemisphere, 
at a position approximately 34° north of the equator. But if this were so, 
the directions of the declinations and dips show that India must have 
rotated through nearly 180° or alternatively have drifted across one 
of the poles to its present position. 

Of these various possibilities the first, that of a northward drift and a 
small clockwise rotation, is evidently by far the most plausible both on 
physical and geological grounds. 

We thus finally conclude that the results are best explained on the 
postulate that India lay about 34° south of the equator during the late 
Secondary and early Tertiary eras, and that at some subsequent time it 
moved into the northern hemisphere and rotated about 25° clockwise 
relative to the earth’s geographical axis. A latitude shift from 34° south 
to the present 20° north, that is, of 54°, corresponds to a linear movement 
of 5000 km. 

It is very hard to see how this conclusion can be avoided except by 
some very special and unlikely hypothesis. A conceivable way of 
reconciling the polarization of these rocks with the present geographical 
position of India would be to assume that the polarizations were all 
originally along the reversed direction of the present field (i.e. south and 
25° upwards, point S in figs. 3 (a) and (b) and that they have been just 
sufficiently unstable to have been all pulled round through some 90° to 
the present position of the peak, that is, 25° east and south and 55° 
downwards. This is so unlikely that it will be provisionally rejected. 

In considering the mechanism whereby such a large movement of 
India might have occurred, one is immediately tempted to consider 
tectonic causes, since current theory of mountain building suggests that 
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India evidently advanced towards the main Asiatic land mass during the 
building of the Himalayas. 

Now the extent of this movement of India, relative to Siberia during the 
Alpine orogeny, is in great doubt and various authorities give figures 
varying from 100 km (Jeffreys 1906) to 3000 km (Van der Gracht 1928). 
Taking a few hundred kilometres as a reasonable mean estimate of the 
relative motion of India and Siberia, the rock magnetic measurements 
would indicate a much bigger northward movement of India relative to 
the geographical axis than relative to Siberia. Whether, in fact, Siberia 
shared in this greater northward movement or not, should be determinable 
by rock magnetic measurements in Siberia. 


§ 5. Five CoMPARISONS WITH OTHER WoRK 


It has been shown by Runcorn (1955) and more recently by Graham 
(1955) that the directions of magnetization found in certain British, 
North American and Icelandic rocks are consistent with the assumption 
that the earth’s mantle has moved as a whole—that is, according to the 
polar wandering hypothesis. 

Of particular significance, by reason of the large number of specimens 
studied, is Graham’s comparison of his results for Permo-Triassic sand- 
stones in North America with the results of Clegg et al. for the Triassic 
sandstones of England. Taken together these results indicate a pole 
position 180 million years ago as 45° N. and 130° E. 

Graham considers that no rotation or translation of Europe relative to 
North America as great as 20° would be consistent with the magnetic 
results. Gold (1955) has given theoretical arguments which indicate that 
polar wandering is much more likely to have occurred than relative motion 
of continents. 

It is interesting to consider our results for the Deccan Trap in relation 
to these hypotheses. The position of the north pole in Kocene-Cretaceous 
times (70 million years ago) as deduced from our measurement is 


28°N. 85° W. 
This can be compared with the estimate by Runcorn of the pole position 


in Eocene-Cretaceous times, as deduced from measurements of British, 
North American and Icelandic rocks, which is 


20 7IN. 1380" He 


It is evident that this large difference between the pole position calcu- 
lated from our results and that calculated by Runcorn indicates that 
India must have moved relatively to Northern Europe and North 


America. ; ; 
Some slight evidence for a rotation of India relative to South Africa 


is obtained by comparing our results with those of Gough (1955) for the 
Pilansburg Dyke system, which is thought to be between 300 and 400 
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million years old. These measurements gave the following values for 
the direction of the polarization : 

Declination N. 24° E. 

Inclination + 69° down. 


The corresponding position of the pole—the north pole if the earth's 
field was reversed or the south pole if it was not—is given by Gough as 
8° north and 43° east, which is situated in Abyssinia. 

Using similar arguments to those which we have applied to our Indian 
results, one can conclude that the most likely explanation of Gough’s 
result is that South Africa has drifted right across the south pole to 
reach its present position and has rotated 24° clockwise in so doing. 
Although the general northward movement of South Africa in the last 
300 to 400 million years is in the same direction as that of India in the last 
70 million years, the relative rotations are in the opposite sense : for South 
Africa 24° clockwise, for India 25° anti-clockwise ; this gives our relative 
angular rotation of nearly 50°. However, the wide differences in ages 
of the rocks makes this comparison of only suggestive importance. 

Whatever the mechanism responsible for the postulated northward 
drift of the Indian peninsular, our present results clearly indicate many 
lines of future research in rock magnetic studies in India. 

In the first place, the discrepancy between the measurements at Linga 
and Khandala require further investigation. It seems possible that some 
part of the difference may be due to the secular variation. To investigate 
this more fully, it is evidently necessary to carry out a more detailed 
survey of the two locations, and to extend the scope of the work to cover 
other parts of the Deccan Trap area. It will be noted that the Linga 
rocks have the smaller dip ; whereas if they are in fact older, they would 
have been expected to have larger dips, due to the postulated northward 
drift and so to their greater southern latitude when formed. 

Secondly, it is well known that Southern India suffered heavy glaciation 
durmg Permo-Carboniferous times. If climatic conditions then prevailing 
over the earth were similar to those of the present day, this would suggest 
that some 250 million years ago India was in a high southern latitude, say, 
not less than 65° south. The latitude of 34° south found in the present 
investigations for Cretaceous-Eocene times is intermediate between this 
and the latitude today. It is therefore important to make further 
measurements on rocks of both pre-Cretaceous and post-Eocene ages 
in order to obtain more detailed information concerning the position of 
India during these periods. This investigation is now proceeding. 
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ABSTRACT 


In this paper p-mesic energy levels are calculated for various nuclear 
charge distributions. Only the Coulomb attraction between the »-meson 
and the nucleus is considered. It is found that the energy levels are 
almost independent of small changes in the charge distribution provided 
the ‘radius of the distribution ’ (defined by the fourth moment of the 
distribution) is kept constant. 

The calculations were done with the aid of a desk machine. The 
methods used are fairly simple and can easily be used with any reasonable 
charge distribution. — , - 


§ 1. INTRODUCTION 


RECENTLY there has been considerable interest in the nuclear charge 
distribution and various experimental data have been interpreted to give 
a value for the radius of this distribution (see, for example, Hill and Ford 
1954b). The most satisfactory of these, both because of the possibility 
of making accurate measurements and because of the comparative 
simplicity of doing calculations, are the differential cross section for 
high energy electron scattering and the energies of x-rays from p-mesic 
atoms. 

In this paper p-mesic energy levels are calculated for various nuclear 
charge distributions. Only the Coulomb attraction between the nucleus 
and the .-meson is considered. Other effects are probably small and are 
discussed by Hill and Ford (1954 a). ; 

Low energy electron scattering (up to ~20Mey) is determined 
completely by the ‘radius of the charge distribution’, R,, defined by 


rn | p(r)rt dr 8ZeR3y, 6.» 3 4 6 ms 0) 
0 
Ri, A2210- om: Ae etm, eat aS 


However, electron scattering experiments have been done more accurately 
at energies from 80 to 180 Mev (Hofstadter e¢ al. 1953) than at low energies. 
From these experiments 7) is found to be in the range 1-16 <r)<1-26. 


t+ Communicated by Professor R. E. Peierls, F.R.S. 
{ Now at Royal Holloway College (University of London). 
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(This can be deduced from the work of Ravenhall and Yennie (1954) ; 
the same result was found by Brown and Elton (1955).) The small and 
medium angle scattering at even higher energies is relatively insensitive 
to the detailed features of the charge distribution. It is interesting to 
see whether further or different information about the nuclear charge 

istribution can be obtained from a study of x-rays from jp-mesic atoms. 

The charge distributions used here are chosen from those used by 
Brenner, Brown and Elton (1954) and by Brown and Elton (1955) for 
high energy electron scattering. They are : 

(i) The square distribution 


p(r)=pp 7<, : 
pir eOOn ere Ray Ry eS Sh 


(ii) Distributions of the form 


1 
P= Po TT exp {= RY eee ee (2) 


with (a) p?=R?/4 and (b) 6?=R?/2. The quantities py and R occurring 
in (3) and (4) are determined from (1) and the condition 


Aa | p(r)r? dr=Ze. OE ee Vf cn pres (0) 
J0 

It is particularly interesting to investigate whether the length 8 of the 
transition range affects the energy levels. In other words an attempt 
is made to answer the question: ‘‘ Does the gross energy-level structure 
of the u-mesic atom give information only about the extent of the charge 
distribution, or can it be used to distinguish between various distributions 
with the same radius R,?” Since the energy of the Is level in lead is 
only about half what it would be for a point nucleus with the same charge, 
it is clear that the difference between the potential of a nucleus of finite 
extent and that of a point nucleus cannot be considered as a small 
perturbation. It might be thought, therefore, that the difference between 
the energy levels corresponding to these types of nuclei would depend on 
higher moments of the charge distribution than the fourth (which determines 
R,, see (1)). However, it will be shown that the energy levels are almost 
independent of small changes in the charge distribution provided &, is 
kept constant. 

While this point was being investigated convenient methods for 
calculating -mesic energies were found. These methods are particularly 
useful for a hand computor. In general it is only necessary to evaluate 
power series expansions for the quantities involved, though for some 
distributions a small amount of numerical integration must be done. 


§ 2. FORMULATION OF THE PROBLEM 


It is assumed that the .-meson is a Dirac particle, i.e. that it has spin > 
(Tiomno and Wheeler (1949), Fitch and Rainwater (1953) believe they 
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observed a 0-2 Mev splitting of the 2p levels for lead, which would confirm 
the assumption), and that its wave function satisfies the Dirac equation} 


{a.ptfu—Vinwn=Bylr) 2 ee ee 8) 
where a, 8 are the usual Dirac matrices, p, » are the momentum and 
mags of the j.-meson. The potential V(r) is determined from the equation 

V72V (r)=4rep(r) 2 2 a oe 
so that V(r)~—Za/r as 7—> 00 
where Z is the atomic number of the atom and « is the fine structure 
constant. Since the charge is concentrated near the origin, V(r) has the 
form —Za/r for all r greater than a certain R which depends on the charge 
distribution but is of the order of the radius of the nucleus. 

Although eqn. (6) can be solved directly, it is of interest to estimate 
relativistic effects. The wave function -4(r) occurring in eqn. (6) is a 
four-component function. It can be replaced by two two-component 
functions ¢,, ¢, satisfying 


(V+uybe+o-Pby 0. gy 
(H—V—p)b,+6-Pba=9, J 
where the o’s are the Pauli spin matrices. The small component ¢, can 
be eliminated to give 
p? B+V p? (W'+V) (H'+V) 
(GetV+#) | 2 Pe ee joe 
1 dV 0¢, Ll V, 
4u2 dr Or © Qu? r dr 
where #’=u—E, and L and S are the orbital and spin angular momenta 
of the ~-meson. The approximation H’+V=—p?/2u may be used in 
the first two terms in the second bracket since they are small. The third 
term in this bracket is much smaller and may be neglected. This gives 


p? p! Lod od. el elle | 
(+0 +8") Popa Pt at Pk Sho. (10) 


L.Sd;=0.. 2s edo) 


The bracketed terms are those of the non-relativistic Schrédinger 
equation, and the last term is the spin-orbit coupling term. 

The kinetic energy p?/2 is roughly 42’. Even in the Is level, where 
it is largest, H’ is only about y/10 (see table 1). The term p*/8u2 thus 
never exceeds about 1:5°%, of the total binding energy. 

To estimate the remaining term it is necessary to evaluate 


Using Poisson’s equation this may be written 


me (® 3Za FP 
[=— YE iE phy? dr~— 32 iF pyr? dr. 


eee 
+ Units such that 4=c=1 are used throughout. 
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R 
Now for the 1s state in lead | bor? dr~t 
0 


(Cooper and Henley 1953), and for other states it is very much smaller. 
Thus this term contributes about 


3Za 
T6u2R3 ~0-2 Mev 
to the ls energy and does not affect the other levels significantly. The 
actual difference between the relativistic and non-relativistic 1s level is 
much less than this (table 1) since this term is about the same size as 
the term p*/8x°, but has the opposite sign. 

Thus the relativistic corrections to the non-relativistic .-mesic energy 
levels are small and, except for the 1s level, the largest contribution comes 
from the term p4/83. This depends mainly on energy and, to the accuracy 
needed here, should not vary significantly for small changes in the charge 
distribution. When spin-orbit coupling can be neglected (s states) the 
energy levels of the u-meson will be given to better than 1% by the 
Schrédinger equation 


{Ve Vehbey BM, ne eee) 1) 


However in cases where it is necessary to include the spin-orbit term it 
is easier to solve the relativistic eqn. (6). 


§ 3. MmerHop oF CaLcuLATING Non-Reativistic ENERGY LEVELS 
It is useful to discuss the non-relativistic problem first as it is somewhat 
simpler in structure and illustrates the essentials of the method of finding 
energy levels. 
The radial wave equation can be written in the form 


2 

Bee te ply PeeCeR Son see (12) 
da ee 

where F/r is the radial wave function and for convenience the following 

change of scale has been made : 


DLO, 
V(r) =— 2 Z?0?¥ (2), eee a ae (13) 
B= 25770702. | 
Rie eae 
Now substituting ago" 
in eqn. (12), 7 satisfies 
724 a see ey meee <= (14) 


In general the function 7 which corresponds to the solution of (12) regular 
at the origin is different from that which corresponds to the solution 
regular at infinity, but for those values of C2 which are eigenvalues of (12) 
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they will be the same. The calculation proceeds as follows: The value 
of the function 7 corresponding to the solution regular at the origin is 
calculated at some point c=NX for various values of C? and the results 
plotted against H(—C?/2uZ?x2). Then the procedure is repeated for the 
solution regular at infinity, and the points of intersection of the two curves. 
give the required p-mesic energy levels. 

In practice it is not necessary to calculate many points on either curve. 
Since the potential V(r) is never greater in magnitude than Za/7, the binding 
energy of a j-meson in any level is never greater than that of a «-meson 
in the corresponding level in the field of a point nucleus. Except in the 
case of the 1s level, where the binding energy is about half what it would 
be for a point nucleus, the difference is not very great and it is less for 
higher levels than for lower ones. It is therefore possible to make a rough 
estimate of the energy levels before starting on a calculation. Suppose 
now that 7, is the function corresponding to the solution of (12) regular 
at the origin for some value C2 of C? and consider the change 67 in 7 
corresponding to the change 3c? in C?. Then, to first order, 57 satisfies 


the equation 
2ngdq-+ 89 =8C. ia, +) eee 


This equation can be integrated to give 


x xX a 
dn(X)=dC2 exp | - } 2” a | | exp | [ 2no(y) iy | dx. (16) 
1 di’, 
Weuae 
where /’) is the solution of (12) Belay: to the value C,? of C?: 
Therefore 


Now No= 


8y(X)=8C? F2dx 


1 : 
(F(X) rR, 
and, provided X is not in the region where F') has become small, é(X) is 
of the same order of magnitude as dC? and is proportional to it. Thus the 
calculation of 7(X) for three values of C? (or #) in the neighbourhood of 
the required energy level provides a check against random errors. in 
calculation since these values should lie on a straight line. This argument 
applies equally to the function corresponding to the solution of (12) 
regular at infinity. The required energy level is therefore given by the 
intersection of two straight lines. 

If the point =X is chosen so that Y (x) has the form We ior 7X 
(corresponding to V(r) having the form —Za/r for r>R=X/2uZa, as 
mentioned in § 2), the solutions of (12) regular at infinity are the Coulomb 


functions \ 4 
F (x)= W 4/20, 141/2(2Cx) 
where W,, »(z) is a Whittaker function. This gives 


a 20\12 W (2C:x) 
n(%)= —C + —— -(=) (41 x) 1/20—1)2, 141 
av 20} Wy), 1412 (202) 
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which can be evaluated using the asymptotic expansions for the Whittaker 
functions when the principal quantum number » is given by n=/-+1, 
nl. In the case n=1, the Is state, 7 can be found directly by solving 
eqn. (15) in a power series in 1/z. When n-4/-+1 the calculation is more 
difficult since in general the asymptotic expansion of the Whittaker 
functions cannot be used and it is then often easier to find the relativistic 
solution. 

The function 7 corresponding to the solution of eqn. (12) regular at 
the origin can be found by solving eqn. (14) directly. When (x) has a 
power series expansion convergent at x=X, it is usually possible to find 
a power series solution of (14) which is also convergent at «=X. 


§ 4. METHOD oF CALCULATING RELATIVISTIC ENERGY LEVELS 


The Dirac eqn. (6) can be separated to give 


1 
{2.—E'— V(r) |F ‘Fz ! “ba=0, 


; ue mee.) 
{—E’—V(r)}@+ \z _ “| F=0 

: dr dest 

where F/r, G/r are the radial components of the wave function and 

E'=p—E is the binding energy and is determined by the condition that 
F and G should both vanish at the origin and at infinity. 

The most convenient units to use differ from those of § 3 by a factor Zx 

and are defined by ae 

Zur, 
E’=2upe, 
V(r)=— Fy). 


Then ¥ (y)+Zz/y as y> oo and the eqns. (17) become 
Liat 

ire ouneee “lano, | 

e+ Vy) F arate ), 


{41 Vy) G4 i “) F=0. 
The required values of «, ie. those for which both # and G are regular 
at the origin and at infinity, can now be found in essentially the same 
way as the eigenvalues of eqn. (12) were found in § 3. Here the ratio F/G 
plays the part of the logarithmic derivative in the non-relativistic case. 
The procedure is to find those values of « for which this ratio for the 
solution regular at the origin and that for the solution regular at infinity 
fit together to form a continuous function. 
Writing F/G=€, eqns. (18) give 
dé 


2 : 
US erect) Schaap war oo eae (LY) 


(18) 
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By an argument similar to that in § 3 it can be shown that for either of 
the required solutions the value of € at any point varies linearly with « 
for small changes in ¢, and a similar method can be used for finding the 
required eigenvalues. 

It is particularly easy to find the 1s level from eqn. (19). The solution 
regular at infinity can be found by substituting a power ser ies in 1/y for €, 
and that regular at the origin by using a power series in y. For the square 
erribarcdl these both converge at the edge of the nucleus, y=2yR,, 
and for the other distributions used the latter series may be continued 
to a point outside the nucleus from a point within its radius of convergence 
by doing a small amount of numerical integration. 

The method used by Gordon (1928) to find relativistic hydrogen-like 
wave functions can be followed through to find the functions / and G 
regular at infinity. It is then found that the corresponding function € is 


e \¥2P(y)+Qy) 
t= (5) "pase OD 
—2y)! 2y) ! 
where PYe)= a bey = or Pky) 


and similarly for Q(y), and 
P(y)=(2by)"F{o—A +1, 20+1; 2By}, 


Qa(y)=— EE (apy Plo, 2041: By} 
y= nr— Ze, 
B?=e(1—e), 
A=Za( vee 
The function F{a, b; z} is a confluent hypergeometric function. Since y 


and A are both ies ot the formula (20) can always be evaluated and 
the difficulties of the non-relativistic case (where the corresponding 
formula contains factorials of negative integers) do not arise. 
The expression (20) for € can be written in the more compact form 
peas. ( € ) W 412, (2BY)—(x+-Ze/2B) W 4_ 1/2, (2Py) (22) 
W 4 +12, (2BY)+ («+ Ze/2B) W 41/2, ,(2By) ” : 
This expression can be evaluated sufficiently accurately using the 
asymptotic expansion for the Whittaker functions W,, ,,(z) for the lowest 
energy level corresponding to each negative except «=—1 (the Is state), 
since then x+Z«/28=0 for the point nucleus and hence «-+Za/28 will 
be small for the values of « used in the calculation. 


l—e 


§ 5. RESULTS 
The methods of the last two sections have been applied to the calculation 
of several -mesic energy levels in lead. The results are given in the tables 
below for a j.-meson of mass equal to 207 electron masses (Smith, Birnbaum 
and Barkas 1953). The potentials used are defined in § 1. 
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The rather strange values for 7) used arise because the calculations were 
in fact performed for a u-meson of mass equal to 210 electron masses and 
the results appropriately modified. 


Table 1. Hnergy Levels of a -Meson in the Field of a Point Nucleus 
and of two Charge Distributions of Finite Extent (see § 1) 


Energies are given in Mev 


Non-relativistic Relativistic 
a O+ Tei 1-217 O+ 1-217 1-217 
B?/R? Of z Of oy 
Is 18-94 | 10-33, | 10-46 21-03 | 10-39 10-53 
2Pse 4-73 4-53 4-84 4-60 4-60 
2Pgols | 14-21 5-81 - 16-19 5-79 5:93 
2s 4:73 5:39 3°55 


Table 2. The Variation of the Non-Relativistic 1s Energy Level 
with Radius for a Square Distribution 


To OF 1-116 1-217 1-319 


ls energy in Mev] 18-94 10-79 10-33, 9-91 


Table 3. The Variation of the Non-Relativistic 1s Level with the 
Length of the ‘ Transition Region ’, 8, for 7)=1-217 


| BR? OF 


rs 
ls energy in Mev 10-33, 10-41 10-46 


7 Point nucleus. { Square distribution. 


tol 


Energy levels for other ».-meson masses can easily be calculated from 
the results given here. Examination of eqns. (12) and (17) shows that, 
since the potential can be written in the form V(r)=f(7/R)/R, with 
R=r,A?10-8 cm, changing the p-meson mass from jt, to Ha Is 
equivalent to multiplying ry by ,/42 and the binding energy by p12//44. 


§ 6. CONCLUSIONS 


Tt has been shown how, given any charge distribution, j.-mesic energy 
levels can be calculated without much difficulty. 

The results obtained show that the 2p—Is transition energy for a square 
distribution of radius 1-21741/10- cm is slightly smaller than the 
experimental value of 6-0 Mev (Fitch and Rainwater 1953). Since the 
1s level is most sensitive to changes in the radius of the charge distribution 
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(it changes about ten times as much as the 2p level for the same change in 
radius), this difference could be corrected by taking a square distribution 
of smaller radius. The 1s binding energy has been calculated for various 
radii (see table 2) and this gives the approximate change of the 2p—Is 
transition energy. From this it may be deduced that a square distribution 
of radius given by an 7, in the range 1-16 to 1-18 would give the 2p—Ils 
transition energy correctly. The calculations of Hill and Ford (1954 a) 
show that it is in fact given by a distribution with ry7=1-17. 

The charge distributions (ii) (a) and (b) defined in eqn. (4) are essentially 
square distributions with a Gaussian ‘tail’. It has been found that 
distributions with varying lengths of tail (the 6 of eqn. (4)), but with the 
same (r?),, give only slightly different values of the ls energy level. 
It is therefore not possible to obtain much information about the nuclear 
charge distribution from energy levels in p-mesic atoms which could not 
be obtained from electron scattering. 
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ABSTRACT 


Earlier work by the Manchester Jungfraujoch group on slow, negative, 
heavy cosmic-ray particles observed in a magnet cloud chamber has 
been repeated. An upper limit to the frequency of slow K~-mesons 
which do not decay in the cloud chamber is now found to be ~1% of 
the frequency of slow protons. From other published data the probable 
ratio of the two frequencies is estimated as about 1: 2000. The identity 
of the observed particles is discussed. The previous results giving a high 
frequency and a long mean lifetime for K~-mesons are not confirmed. 

One event showing the possible simultaneous production of two slow 
K--mesons and one example of an upward-moving S*-particle are 
described. 


§ 1. INTRODUCTION 


RECENT work with cloud chambers and nuclear emulsions has led to the 
identification of several species of positive K-meson. The identity of the 
negative K-mesons, on the other hand, is still uncertain. This is not 
only due to the relative infrequency of the K™-particles but also to the 
methods of observation. In emulsions and multi-plate cloud chambers 
both the nature and the sign of a decaying particle can only be determined 
unambiguously if the charged secondary or secondaries are 7-mesons 
which can be followed to the end of their range. So far, no decay events 
of this type have been established for negative particles, most of which 
interact with atomic nuclei at the end of their range (cK events). In 
such events, however, no information on the nature of the K-particle 
ean be obtained. 

The results from magnet cloud chambers are therefore of particular 
interest. In such work only the 7~-meson has hitherto been identified 
with certainty (van Lint and Trilling 1953). There is also evidence for 
the y~-meson (Armenteros e¢ al. 1955 a) and in addition the observations 


+ Now at Imperial College, London. 
+ Communicated by Professor G. D. Rochester. 
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of Buchanan etal. (1954) indicate the existence of negative particles 
decaying with a mean lifetime exceeding 10~* sec. Most of the results 
of cloud chamber and emulsion work suggest a low frequency for 
K--mesons. However, the experiment of Astbury et al. (1953) apparently 
showed the existence of a negative K-particle occurring with a high 
frequency and having a mean lifetime exceeding 4x10-*sec. The 
frequency of slow ,K--mesons relative to protons with approximately 
the same velocity found by Astbury et al. exceeds the frequency of all 
K-mesons found by Dahanayake et al. (1955) in a systematic study of 
‘grey’ tracks from nuclear interactions of energy greater than ~2Gey. 

Since the results of Astbury et al., if confirmed, would also suggest 
that not all K~-particles give rise to oK-events when brought to rest, 
we decided to repeat the experiment in order to obtain a better estimate 
of the frequency and, if possible, of the lifetime of the particles observed 
by them. 

Our method consisted of finding tracks due to heavily-ionizing particles 
whose origin could be located and whose direction of flight and sign 
could therefore be determined. Tracks indicating particles with negative 
charge were further examined. Those due to light mesons were 
eliminated by momentum measurements and visual estimates of 
ionization. By choosing tracks of sufficient length the probability 
of including negative hyperons was made small, and it was assumed 
that the remaining tracks were due to K~-mesons. By correcting for 
the inefficiency in scanning the photographs and for the expected 
number of spurious events due to upward-moving protons the frequency of 
slow K~-particles was estimated and compared with that of slow protons. 

We attempted to estimate limits to the frequency of the K -particles 
by considering statistical fluctuations in their rate of occurrence and by 
allowing for extreme values of the scanning efficiency and of the number 
of spurious events. We also considered the nature of the particles 

observed by us and the possibility of estimating their mean lifetime. 
_ No attempt was made to assess the frequency of K*-mesons since the 
uncertainty in visual estimates of ionization prevents the separation of 
a small number of K*-mesons from the much larger number of protons. 


§ 2. EXPERIMENTAL DETAILS AND RESULTS 


The magnet cloud chamber at the Jungfraujoch has been described 
by Astbury ef al. (1952) and by Newth (1954). When the photographs 
considered in this paper were taken, the chamber contained a wedge- 
shaped lead block immediately above the illuminated volume and a 
12mm copper plate across its centre. Since the earlier investigation of 
Astbury et al. (1953) the depth of the chamber has been increased so as 
to double its illuminated volume. The photographs used by us were 
therefore better suited to the observation of particles coming from well- 
defined origins. The results presented here were taken from 13 000 
photographs, all of which were scanned at least twice. 
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The following criteria were used in the selection of possible K~-events: 

(1) The track must show evidence of a likely origin. 

(2) Its curvature must be uniform and indicate a negative charge 
assuming that the particle travelled downwards from its supposed 
origin. 

(3) The ionization density of the track must be at least twice that 
of the track of a fast particle (ie. [>2.x J,). 


(4) The track must be more than 10 cm long. 


(5) Neither the photograph containing the track nor any neighbouring 
photograph must show evidence of any gross convective distortion 
occurring in the cloud chamber. 


In establishing a likely origin it was demanded that, on three- 
dimensional reconstruction, the track should pass within 0-5cem of 
a nuclear interaction. The limit of 0-5cm was based on experience 
with the stereoscopic reprojector. The stipulation of an origin allowed 
the rare, downward-moving, negative particles to be distinguished from 
the comparatively frequent, upward-moving protons. 

The requirement of uniform curvature was made to exclude tracks 
whose curvature is falsified by distortion at the edges of the illuminated 


Measurements of Tracks of Slow Negative Particles 


Serial Particle Estimated Mass of | 


number momentum ionization particle 
of event (mMev/c) (Io) (Mme) 
SJ 672 ans 3 to 6 1050 to 2700 
SJ. 890 170425 >10 >1100 
SK 560 2807 8° 2to4 500 to 1600 
SQ 1437 A Loo 4 =i >850 
SQ 1437 B Oe, >6 >700 


ee ee 


Measurements made on five tracks apparently due to negative heavy particles. 
The second column gives the momenta of the particles, the third the ionization 
density in the track and the fourth gives the mass values obtained by combining 
the figures from the previous columns. Event SQ 1437 is reproduced in PI. 19. 
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volume of the chamber. The uniformity of the curvature was established 
by microscope measurements of the track. Tracks at least 10 cm long 
were demanded to ensure an accurate measurement of their curvature 
and to reduce the likelihood of their being due to negative hyperons. 
The mean lifetime of charged hyperons is less than 10~1°sec (Davies 
et al. 1955) and the mean decay length for heavily-ionizing hyperons is 
therefore less than 1-5 cm. 

An ‘origin’ was defined as a nuclear interaction satisfying the 
following criteria : 


(1) It must be situated in the lead wedge or copper plate and must 
lie within 6 cm of the illuminated volume of the chamber. 


(2) In addition to the negative particle it must give rise to either 
ewo or more fast particles other than electrons or at least one heavily- 
ionizing particle. 


Four events were found which fitted the criteria listed. In no case 
did the slow negative particle decay in flight in the cloud chamber. 
Details of the events are given in the table. The momenta quoted 
have been obtained and corrected by the method described by Barker 
(1954). The ionization estimates are the combined results of visual 
estimates by six observers ; they are corrected for the dip of the tracks 
where necessary. 


§ 3. Discussion oF EVENTS 


The selection criteria listed in §2 do not take into account the 
possibility that an upward-going proton without a visible origin may 
accidentally pass so close to a nuclear interaction in the plate or wedge 
that it appears as a negative particle coming from it. To estimate the 
number of such spurious events expected in 13000 pictures we deter- 
mined the number of tracks of upward-moving slow protons which 
showed no origin but otherwise fitted our selection criteria. From this 
number and the number of nuclear interactions observed -we conclude 
that the expected number of spurious events in the sample is ~1 for 
the lead wedge and ~0-1 for the copper plate. 

It follows that the probability of both events from the plate (SJ 672 
and SK 560) being spurious is ~0-01 and we therefore consider that 
at least one of them is genuine. Event SJ 890 from the wedge may 
be spurious in view of the greater probability of a chance coincidence 
in the wedge and of the high estimate obtained for the mass of the 
particle. 

Event SQ 1437 (PL. 19) is of great interest. The two steeply dipping 
tracks marked appear to come from a common point in the wedge. 
Mass estimates show that they may be due to negative K-particles. | 
The point of intersection of the tracks is defined to within -0-2 em. 
Since it lies only 1-5 cm in front of the plane marking the rear surface of 
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the illuminated region of the chamber it is possible that associated 
particles, if any, would not have been observed. From the measured 
rate of upward-going protons entering the wedge we estimate the proba- 
bility of a chance association of two such particles to be ~5 x 10-3 for a 
sample of 13 000 photographs. 

Two events tentatively identified as the production of K*—K~ pairs 
have so far been reported (Davis ef al. 1955 and Ceccarelli et al. 1955 b). 
As interpreted these events are consistent with the formalism proposed 
by Gell Mann (1955). The same formalism would require the production 
of at least two other unstable particles to explain our event in terms of 
two K”-mesons. However, the possibility that our event may be spurious 
cannot be ruled out. 


§ 4. THE RELATIVE FREQUENCY OF K~-Mrsons AND PROTONS 
(i) Hstimate of Frequency 


For comparison with the results of other workers we have estimated 
the ratio of the frequency of our K~-particles to that of protons which 
have approximately the same velocity and fit our selection criteria. 
For this purpose it is sufficient to assume that all slow, positive, heavy 
particles are protons. From the number of protons found in a sample 
of 1000 pictures the number in 13000 pictures was calculated to be 
1070+130 from the copper plate’ and 1900-+180 from the lead wedge. 
In this sample the scanning efficiency for slow protons was high and we 
assume it to have. been unity. We estimate the efficiency of detecting 
non-decaying K~-particles in the batch of 13000 pictures to have been 
about 0-5. 

Assuming both events from the plate to be genuine we obtain a value 
for the frequency of slow K~-particles relative to that of slow protons 


N(K-)/N(P)=(2/1070) x (1/0-5)~1/250. 


The velocity ranges of the K-particles and the protons considered by us 
are slightly different (0-3<8<0-6 and 0-1<f<0-6 respectively). 

To estimate the upper limit of the frequency, the lower limit to the 
scanning efficiency is taken as 0-3 and all events which fit our criteria 
are considered as genuine. The possibility that the number of events 
present represents a fluctuation from the true mean rate is taken into 
account at the 5°% confidence level. Values for the upper limit may be 
obtained for events from both the wedge and the plate. They are 


from the lead wedge N(K7~)/N(P)<1/100, 
from the copper plate N(K~)/N(P)<1/80. 


No useful lower limit may be set since the probability of both events 
from the plate being false is ~0-0I. 
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(ii) Discussion 

This method assumes that the angular distributions of protons and 
K--mesons are similar. It is evident that the results obtained from 
experiments with nuclear emulsion are free from this possible source 
of bias. 

The upper limit to the frequency of K~-mesons disagrees with the 
value obtained by Astbury e¢ al. (1953)} but is in agreement with other 
work. Armenteros et al. (1952), using a magnet cloud chamber, found 
one non-decaying, slow K~-meson among 900 protons. Since the 900 
protons came from production layers both within and above the cloud 
chamber their result is not incompatible with ours. 

In photographic emulsion Dahanayake et al. (1955) found ten decaying 
K-mesons among 2320 protons by tracing all grey tracks from cosmic 
ray stars to the end of their range. The results of Hornbostel and Salant 
(1955) and Chupp et al. (1955) indicate that at least 90°, of artificially 
produced K~-mesons interact when coming to rest. If the K”-mesons in 
the cosmic radiation behave similarly the ten decaying K-mesons found 
by Dahanayake et ai. would all be positive. Ceccarelli et a/. (1955) in a 
systematic search for slow K-mesons found a value of 10-3 for the 
ratio of the number of positive and negative particles. It follows that 

N(K-)/N(P)=N(K-)/N(K+) x N(K+)/N(P)~1/10 x 10/2320~1/2000 
where the numbers refer to slow particles observed in nuclear emulsion. 

The observed N(K)/N(P) ratio will not be representative of the ratio 
of frequencies in other energy ranges unless the energy spectra of the 
two types of particle are similar. There is good evidence that the 6° and 
proton spectra differ at energies below about 200 Mev (James 1955), but 
it is likely that the production mechanism of K~-mesons is different from 
_ that of the majority of 6°s (Gell Mann 1955). We therefore cannot 
assume the K~-spectrum to resemble either the proton or the 6° spectrum. 

To obtain an accurate measurement of the frequency of K~-particles 
by our method a very large number of photographs would have to be 
examined. Such a search would not be worth while unless the systematic 
errors due to uncertainty in the scanning efficiency, spurious events, etc. 
could be reduced. 


+ We have re-examined the 10 K--tracks found by Astbury ef al. and find 
that only two fit our stricter criteria. Moreover, it has been pointed out to us 
by the authors of this work that there is an error in the calculation of the 
number of spurious events due to upward-moving protons. This arose from’ 
using the average number of nuclear interactions observed on the photographs 
to calculate the probability of a chance coincidence in space between a proton 
track and a nuclear interaction. If the criteria of association are not very 
strict it is possible for a single photograph to show so many nuclear interactions 
and so many proton tracks that a chance association is almost certain to arise. 
To estimate accurately the number of these chance associations it is necessary 
to consider the frequency distribution of nuclear interactions and proton tracks 
on the photographs. We have made a careful check of this point in our work 
and have found that, with our strict criterion of association, it is adequate to 
use the average numbers quoted above. 
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§ 5. Ipentrry or tHE K~-Parricies 


The following modes of decay have been proposed for positive charged 
K-mesons : . 


Ka pte 
K, p+, 
Kp oer ue. 


T>a+r+n, 
7 >7+7°+79, 
X>7+7°. 


The number of Kg-decays observed is small compared with the numbers 
of other K-decays present in photographic emulsion and none has yet 
been reported in cloud chamber work. Armenteros et al. (1955 b) have 
shown that the negative K,-meson is either very rare or non-existent, 
but in view of the difficulty of identification we think it possible that the 
K™, may yet be frequent in comparison to other K~-particles. 

Observations of decays at rest in stripped emulsions suggest either 
that the «-meson is rare or that it has a very short lifetime (see Alvial 
et al. 1955). 

A single magnet cloud chamber is inefficient for observing decays of 
long mean lifetime such as that of the 7-meson. Hence the low frequency 
with which 7-decays in flight are observed (9 in 100 000 photographs by 
our group) does not mean that non-decaying r-mesons are rare. The 
existence of the negative 7-meson is well established and some or all of 
our events could therefore be so identified. 

The situation with regard to y -particles is somewhat less clear, 
but direct evidence for their existence has recently been obtained by 
Armenteros .et al. (1955a) and Arnold et al. (1955a). There is, in 
addition, indirect evidence from V -events, some of which can most 
easily be identified as y -particles (Trilling and Leighton 1955). On the 

other hand, Alvial et al. (1955) have shown that slow y*-particles are 
more abundant than slow r*-mesons and Bridge et al. (1955) have 
found the mean lifetime of y‘-mesons to be greater than 10~° sec. 
Therefore the identification of our events as y -mesons is also possible. 

In view of the scarcity of evidence on the nature of K--particles with 
the exception of the easily identified +-mesons there is no strong reason to 
assume that the relative frequencies of the different modes of decay of 
K--mesons are different from those of the K*-mesons. Hence it is 
difficult to draw any definite conclusions on the identity of our non- 

' decaying negative particles. 


§ 6. Discussion oF LIFETIME 


In view of the uncertainty as to the nature of the particles observed 
by us we are not justified in estimating a lifetime by the method of 
Astbury et al. We also point out that even by applying this method we 
obtain for the lower limit to the lifetime a value that is no greater than 
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10-®sec. Buchanan et al. (1954) found a similar limit by the direct 
study of negative charged V-events. An estimate of the mean lifetime 
of artificially produced negative K-mesons based on decays in flight is 
given by Goldhaber et al. (1955) as (0-65-£0-45) x 10~* sec, but Arnold 
et al. (1955 b) find evidence for a short-lived component among )~-events 
not due to hyperons. 


§ 7. A PROBABLE S-EVENT 


In the course of our work we found one event which does not fit all 
our selection criteria but which is apparently an upward-moving, positive 
K-meson which comes to rest and decays in the copper plate. The event 
is reproduced in Pl. 20 and the measurements are given in the caption to 
the plate. While we recognize that the event may be a small nuclear 
interaction in the plate in which a negative K-meson is produced, we do 
not regard this interpretation as likely and we note that the number of 
upward-moving protons stopping in the plate is ~200. We should there- 
fore expect to find about one stopping K+-meson. 


§ 8. CONCLUSIONS 

An upper limit of ~1°% has been found for the ratio of the frequencies 
of K~-particles with velocities 0-3<$<0-6 to protons with 0-1<f<0-6. 
This limit is compatible with the results of other groups. A most probable 
ratio of ~1/2000 has been deduced from the available evidence. No 
definite conclusion as to the identity of the particles has been drawn. 

The previous results of Astbury et al. (1953) concerning the frequency 
and mean lifetime of these particles have not been confirmed. 

One event has been found which can be interpreted as the simultaneous 
production of two K~-mesons, another is probably an S-event due to an 
upward-moving, positive K-meson. 
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ABSTRACT 


A complete treatment is given of the steady state spiral both in growth 
and evaporation of crystals, taking account of the strain energy around 
the dislocation. It is shown that when the crystal is evaporating there 
is a limiting value of the undersaturation beyond which there is no steady 
state solution of the problem. The core of the dislocation is then unable 
to remain closed and a macroscopic pit appears at the termination of the 
screw dislocation. ‘ 


§ 1. IvrropucTION 


ALTHOUGH a theory of the spiral on growing crystals was developed a 
number of years ago by Frank (1949), and Burton, Cabrera and Frank 
(1951), there has been no attempt to extend the treatment to evaporating 
crystals. (The word evaporation is here understood in its general sense, 
that is to say, including dissolution, chemical etching, oxidation, etc.) 
Such an extension is particularly important in view of its possible applica- 
tion to the understanding of the formation of etch pits at dislocations. 
This paper develops such a theory and gives as well a more complete 
treatment of the steady state spiral in growing crystals. To distinguish 
between growth and evaporation it is necessary to take account of the 
supplementary energy localized around the dislocation itself. As a first 
approximation we consider only the strain energy, although other factors 
such as the effect of impurities must be taken in account if we want to 
understand completely the formation of etch pits. We limit ourselves 
of course, to dislocations finishing at the surface with a screw component, 
and giving rise to steps on the surface; the extension to pure edge 
dislocations has to be made from a different point of view and will be 
considered in a separate paper. 

In §2 a fundamental expression for the velocity of advance of a step 
on the crystal surface will be introduced and a discussion will be given 
of the conditions under which it is expected to be correct. In $3, the 
theory of the steady state spiral will be developed both for growing and 
evaporating crystals. Finally a discussion of the physical implications 
of the results obtained will be given in § 4. 


Eee 
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§ 2. THE MovEMENT or STEPS 


The basis of the treatment of the spiral given in $3 is the following 
expression for the normal velocity of advance v of an element of step 
on the crystal face, located at a distance r from the point of emergency 
of the dislocation and having a local radius of curvature p : 


Deel (Pal) Appell tal ns we. (1) 
The quantities v,., p, and ry have the following meaning : 


(1) v is the velocity of a straight step far from the origin. 


(ii) p, is the radius of the critical nucleus for two-dimensional nucleation 
on the crystal surface, given by 


Pe VS) Zee, ce fe, rk (2) 
where y is the surface energy per unit surface of the edge of the step, 
2 the volume per molecule in the crystal and Ay, the difference in 
chemical potential of a molecule in the surrounding medium and in the 
crystal (supposed perfect). According to whether dy, 2 0 the crystal is 
growing or evaporating. 

(iii) Finally rp is given by 

: fj =ub*/Sa*y, So nth, Tan eR deta ee ee 
where pz is the shear modulus and b the Burgers vector of the dislocation. 
79 represents a measure of the strain energy around the dislocation ; 
it is assumed that the dislocation is normal to the crystal surface. 

A justification for the right-hand member of (1), can be given in the 
following manner: If the advance of the step is controlled by some 
diffusion field, the velocity of advance should be proportional to the 
difference in chemical potential in the surrounding medium far and near 
the step. Assuming equilibrium to be maintained near the step this 
difference is Ay, if the step is straight and far from the origin ; if this is 
not the case, there are two main contributions which will change this 
difference by altering the chemical potential near the step: (i) the 
curvature 1/p of the step increases the latter by y2/p=(p,/p) 4, (ii) the 
presence of a dislocation normal to the surface of the crystal introduces a 
density of strain energy roughly equal to pb?/87?r? and consequently an 
increase of the chemical potential b?.Q/87?r?=(rop,/7T?) Muy. These two 
effects decrease Ju, by the factor appearing at the right-hand side of (1). 

It will be noticed that (1) is valid for both growth (4) >0) and evapora- 
tion (4y<0). In the latter case the radii of curvature p, and p become 
negative, this does not change the sign of the second term in (1) but 
changes the sign of the third term because the strain energy due to the 
presence of the dislocation opposes growth and favours evaporation. 
We must also mention that, in taking account of the strain energy during 
a steady state process, we must consider the strain at the interior of the 
erystal, as the strain near the surface will continuously reproduce itself. 
Finally the length, 7), has no physical meaning for small b and large y 
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being then smaller than the atomic dimensions ; however, if b is large 
and y small, Frank (1951) showed that the dislocation should have a 
hollow core of radius given by (3) where now y is the surface energy of 
the surface of the hollow, which is not necessarily equal to our y. In 
what follows the two y’s will be assumed equal, which is a good simplifying 
assumption for close packed surfaces. 

Finally, let us emphasize that eqn. (1) is valid only under the following 
conditions: (i) the step is far away from any other step so that the 
diffusion fields are independent of each other ; (ii) the curvature of the 
step is everywhere small, so that it cannot advance by surface nucleation ; 
(iii) the influence of the crystallographic orientation of the step is com- 
pletely neglected. 

Conditions (i) and (ii) are satisfied only for low | 4uo| (large | p,|)- 
Condition (iii) should be satisfied provided equilibrium around the step is 
maintained ; this is likely to be true for atomic steps advancing under any 
kind of diffusion field. There are also cases, where the condition seems to 
be satisfied for a multi-atomic step centred on a dislocation of large 
Burgers vector (for instance the spirals in SiC crystals); our treatment 
will include these cases. The reasons for these multiatomic steps to 
remain together instead of separating into many atomic steps will 
not be discussed here. . 


§ 3. THEORY OF THE SPIRAL 


Going back now to our problem, it is well known that a step should 
wind itself in a spiral around the point of emergence of the dislocation, 
and that the steady state of this spiral should be represented by a rotating 
one with the equation 


= f(r) cb; 2k, ee ee ne 


in polar coordinates, where w is the angular velocity of (let us say) 
clockwise rotation of the spiral. In terms of w the rate of growth (or 
evaporation) # of the crystal and the distance d between successive turns 
of the spiral far from the centre are given by 


R=b'/27, d=2rv,fm, . 2. . 5 2 5 


where 6’ is the height of the step as well as the component of the 
Burgers vector of the dislocation normal to the surface. 

Our problem is then to find the function f(r) and the value of w, by 
using the expression (1). In polar coordinates, the normal velocity v of 
advance of the step and its radius of curvature p are given by 

v= —rb/[1+ (r6')2]¥2, p= [1+ (r6’)29/2/[26’ +70" +7293], 
where ; 
6=00/0t=—w, 6'=00/dr. 
Replacing v and p in (1) we obtain the equation 


rol L + (70’)?]=0.0[1—(rop_/?2) 1+ (102 B2—v,,p,[20’ +10" +726'3], 


a 


Dislocation Theory of Evaporation of Crystals 453 
By introducing the new non-dimensional variables s, f given by 
t=Sp,, 1(00/er)=s(00/ds)=tand, . . . . . (6) 
where ¢ is clearly the angle between the radius vector and the tangent 
to the step, the differential equation reduces to 
$'=[1—(s,/s?)] sec 6—(1/s) tand—w,s, . 2. . . . (7) 
where the notations 
1o=SoPo> W=WVoo/ Pes SS er ele EAD & ted) 
have been used. The non-linear eqn. (7) has to be solved for d(s) and 


w, given the boundary conditions at both large and small s. In terms 
of w,, R and d are given by 


R>a,b'v,,/27p,, d=(27/w,)p,. — .. . ». . (9) 
In all the cases we are going to consider, it will be assumed that the 


crystal is so large that the step becomes circular at large s ; developed 
in inverse powers of s the solution is easily shown to be 


(s)=$7— 1/(w1s)+1/(w18?)—..., s>1. . . . (10) 
At small s we must consider several cases according to whether s, is 


negligible or not, and in the later case, whether growth or evaporation 
is occurring. 


3.1. 8)=0: Growth and Evaporation 

For small Burgers vector b and large y, the length ry is of the order of 
atomic dimensions, therefore s, is negligible particularly for the large 
values of p, (10-+.cm or more), in which we are specially interested. We 
begin then by assuming s,=0 altogether, which implies that there is no 
difference between growth and evaporation. 

Under those conditions, the step should reach the centre s=0. The 
solution for s small is then easily found to be 


O(S\=258— 5048" 15. ti SSL, ys SD) 
giving the right curvature 1/p, at the centre. The corresponding value 
a),=0-33 is then obtained by matching the two solutions (10) and (11). 
The complete solution 4(s) is represented in fig. 1 (sj=0) and the corres- 
ponding spiral is shown in fig. 2. It should be noticed that the distance 
between successive turns reaches very rapidly the limiting value d=19p,. 
The radius of curvature of the step becomes very rapidly large with 
respect to p, so that for purposes of calculating; for instance, diffusion 
fields the spiral can be replaced by concentric circles distant d from each 
other. 

This solution assumes however, that the point of emergence of the 
dislocation remains anchored as the spiral turns around. On the other 
hand, the step finishing at the dislocation has an energy per unit length 
(tension t=yb’) which is not small with respect to the tension 7’ of the 
dislocation itself (see Nabarro 1952); consequently it should tend to 
make the dislocation move. When growth occurs near the melting 
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Representing the angle 4(s) between the radius vector s and the tangent to the 


parameter 8, (see text). 
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point it might be assumed that the dislocation can climb as well as glide ; 
then, assuming that the dislocation is anchored at a distance h beneath 
the surface (say 4~10~4cm) the point where the dislocation meets 
the surface should turn in a circle of radius r’ of the same order as h, 
or more precisely r’=ht/2T[1—(t/27)?]12 if image forces are taken in 
account. The end of the step should remain tangential to the radius 
vector. Therefore if we define r’=s'p, the solution near s’ should be 

| $(s)=H(s—s’)—hurg(s—s' P+... 

and the value of w, would change with s’. For a value s’=1 we find 
w,=0-28, which does not differ very much from the value for s’=0. 
If growth occurs at lower temperatures, at which only glide can occur, 


the dislocation should move in its gliding plan and the point of emergency 
should oscillate around s=0. 


3.2. 89 >0: Growth 

For large Burgers vector b and/or small y, sp is not negligible. Taking 
the characteristic values y~5 x 10? erg/em?, w~5 x 10!" erg/em*, and a 
small rate of growth corresponding to duy~10-' erg, s, is larger than 
0-1 if 6 is larger than 10~* cm which is known to occur in polytypes. 
On the other hand, the dislocation should have at equilibrium a hollow 
core of radius 7) >10~-° em, which would make it immobile. 

For a given value of s,>0, we want a solution of (7) representing a step 
which becomes circular as it approaches a limiting radius c (as yet 
undetermined and not equal to s,) representing the radius of the hollow 
core of the dislocation during growth. The boundary condition for 
s>c must be then ¢ > —4z, the minus sign because the step is now 
turning in the reverse sense to what it does at large s. For every value 


of w,, there is a solution of (7) which cuts 6=—47 at a certain value 
c of s. All those solutions however, are of the form ¢=—4z7-+a(s—c)1/? 
unless ¢ is the positive root of 
Co eg OM etal a tec «nl Bi hk 2) 
In this case the solution is of the form 
d=—tr+a(s—c), soe, . . .1. . . (13) 
with a given by ca?+ w,ca+ 1+ 2c=0. 


This latter solution is the one we are interested in, as it is the only one to 
satisfy the further condition 6+ — oo as s-c, necessary for a hole to 
exist. The radius of this hole r’=cp, has the following value: 

To =the +4(ro/pe)?—-1y;  - - - - - (14) 
it tends to ry for very large p, and decreases slowly as p, decreases. 

The complete solution ¢(s) for any value of sy, and the corresponding 
value &,(s)), can then be found by numerical integration. Figure | 
includes the solutions for s;=4 and for the rather extreme value s)=1. 
Figure 3 gives the curve (8 ) (8)>0) which is seen to decrease very 
slowly as Sp increases. 
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Angular frequency w, of rotation of the spiral in units v../p,, for different values 
of the parameter s,) (see text). s,<0 corresponds to evaporation, 
8) >>0 to growth. 


3.3. Sy <0: Evaporation 
It is then useful to redefine s and sy as positive quantities, that is to 
say r=|p,|s. Equation (7) becomes 


fh’ =[1-+ (s9/s?)] see 6—(1/s) tan d—w s. 
The boundary condition for sc (¢ is again still undetermined) is now 
. 5 
¢—> tz as the step does not change its sense of rotation with respect to 
s large. The solution for sc is then given by 
p=tn7—a(s—c), sc, ae Soe aes ee oer RL 
provided the value of ¢ is determined from the equation 
6? e-sn0, ee, eT 
and a is calculated from 
122 Cr ce 
c*a*—w,c°a-+ 2c—1=0. 


The two roots of (16) are now positive (provided sy<+, see below) ; 
however, it will be shown in §4 that only the smallest root corresponds 
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to a minimum for the free energy of the evaporating crystal. The radius 
19 =C| p,| of the hollow core of the dislocation becomes then 


Peera Petes | Pel I2), (17) 
which is, of course, the same as (14) if p,=—|p,|. The important point 
in this expression is that for | p,|<47, there is no real solution of (16), 
and therefore no steady state spiral. The physical implications of this 
result will be discussed in § 4. 

In fig. 1, several solutions (s;<0) are given, including the limiting 
one for ss=}. The curve w,(s9) is represented in fig. 3 (s)<0), it is seen 
that w, increases very rapidly, as sy, >}. At sy=1, w,(4)=1-65 and 
joins, with a vertical tangent, the curve w,’(s)) for the solutions corres- 
ponding to the second root of eqn. (16). 


§ 4. Discussion 


The physical reason for the formulae (14) or (17) giving the radius of 
the hollow core of the dislocation does not appear clearly in the argument 
of §3. The following, more general derivation (Cabrera et al. 1954) 
might clarify this point. Let a crystal with a hollow dislocation of 
radius 7 be growing from or evaporating into a medium, the difference 
in chemical potential of a molecule in the medium and in the crystal 
being Ju). An increase of r to r+dr will produce a change dG in the 
free energy of the crystal per unit length of dislocation given by 


AG = (2rr/ 2) Ap, dr+ 2ry dr—2rr(ub*/87r?) dr 
= (2ryp,/r)(c?--e—89) ar, 


where the notations used in §3, r=cp,, 7)=8 p, have been introduced, 
and the explicit assumption that y remains constant has been made. 
The roots of (12) or (16) make, therefore, G stationary. It is also seen 
that the positive root of (12) and the smallest root of (16) make G a 
minimum. It is therefore concluded that the radius ry’ of the hollow 
core of the dislocation remains at the value (14) or (17) while the crystal 
grows or evaporates. - 

The more interesting result in this paper is the existence of a critical 
value | 4u.| for | 4uo| beyond which a dislocation in an evaporating 
crystal cannot remain closed ; otherwise, it must open-up ’. Then the 
crystal evaporates from the free surface as well as sidewise from the core 
of the dislocation outwards, forming a macroscopic pit. At the critical 
value itself (| p,|=4r,) the radius of the core r9’=2r) remains submicro- 
scopic. The critical value | Apy.| is given by 


Appa 2atyeiubemwe x ee 2a (18) 


where (2) and (3) have been used. In the simplest cases of evaporation 
into a vapour or dissolution into a diluted solution we can write 


d = 
SER. 8, VOL. I, NO. 5.—MAY I950 2H 
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| tig |=aer In (po/p) where po is the saturation pressure (concentration) 
and p the actual one occurring above the crystal ; then we obtain 
P-|Po=eXP [— 2777? Q/kT yb? ], 

where p, is the critical pressure below which the dislocation should open-up. 
To illustrate the numerical values that we should expect let us assume 
y~5 xX 10? erg/em?, w~5 x10!" erg/em?, Q~10- cm3, kT’~10~* erg. 
Then we obtain (py—p,)/Porl0? if b~10-® cm and (p,/po)~10* if 
b~10-8 em. This indicates that a very low undersaturation will be 
sufficient to open-up a dislocation of large Burgers vector; on the 
contrary very large undersaturations might be necessary to open up 
dislocations with an elementary Burgers vector. The actual values of 
P,/ Po Will of course be very sensitive to the assumed value for y. 

Another result worth further discussion is the value obtained for the 
semi-vertical angle 8 of the conical surface of the crystal. As long as 
[4.9] is small, it is well known that the surface remains practically flat 
so that f is practically equal to $x. This is also true for the obtainable 
values of du, above a growing crystal. In the case of an evaporating 
crystal it is interesting to consider the value of 8 when the critical value 
| Muy,| is attained. Then 

tan B=d/b=0-19 x pb/y. 
Assuming the same characteristic values for 1 and y we get tan B~2 x 108 6 
(> being measured in cm) which indicates that 6 remains very nearly 
37 for a large Burgers vector but diminishes to values of the order of 
7/3 when b~10-* cm. 

The application of the results of this paper to dislocations of small 
Burgers vector, when the critical value | 4jo,| is approached, should be 
considered doubtful. Then the value of p, is so small that the basic 
expression (1) should not be valid anymore, particularly because the 
internal part of the spiral should then be able to advance by surface 
nucleation. The importance of surface nucleation in the formation of 
etch pits is also implied in the fact that etch pits do form also at pure edge 
dislocations. This point of view will be developed in a forthcoming 
paper. 
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XLVI. Cyclotron Resonance under Anomalous Skin Effect Conditions 
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SUMMARY 

The purpose of this note is to discuss the results to be expected from 
attempts to observe cyclotron resonance in metals. A brief summary is 
first given of the simple theory for a free-electron model, assuming the 
classical skin-effect equations to hold. Theoretical results are then given 
for a free-electron model under anomalous skin effect conditions, and a 
qualitative discussion is given of the results to be expected for a real 
metal under these conditions. 


THE theory of cyclotron resonance is most simply formulated if the r.f. 
fields within the conductor are assumed to be circularly polarized, so that 
e.g., H=H,+iH,=E, e*”’ (where H, may be complex), and the steady 
field H is applied in the z direction. Then for a degenerate gas of quasi- 
free electrons of mass m* and relaxation time 7, it is easily shown that the 
conductivity co, is given by 
Fg=0,/(1+iw’s) eee eel eee (1) 
where o, is the d.c. conductivity Ne?7/m* and w’=w—wo=w—eH/m*c. 
This expression holds under ‘ classical’ conditions, when the mean free 
path of the electron is small compared with the skin depth. If, instead 
of being circularly polarized, the applied electric field is linearly 
polarized, so that H=H,—R(E,e**') the resultant currents are given by 
J ,=R(o,,e""), J, =R(opHh,e*"), where 
op =310¢(w¢e) +o¢(—@¢)], o7= — 3115 (Me) —%(—e)], - ++ () 
and the values thus found for o,, and o7 agree with those given by Dressel- 
haus, Kip and Kittel (1955 a,b). If on the other hand the current is 
linearly polarized, J=J,—R(J,e), then L,,—R(p;J,e""’), etc., and 
Pr=4lpc(@e) + Pc(—@c)]. Pr=— 241 Pele) —Pcl—e)] 5. (Po= I/og), (3) 
and using (1) and (3) we have 
Pr poll -iwt), pp pooner en - a (4) 
If now a specimen obeying these ‘ classical ’ equations is of thickness 
large compared with the skin depth the quantity observable experi- 
mentally is the surface impedance Z=H(0)//,. where (0) is the electric 
+ Communicated by the Author. 


2H2 


460 R. G. Chambers on Cyclotron Resonance 


field at the surface and J, is the total current flowing in the skm. We 
consider three simple cases : 


(a) If the specimen is placed with its surface in the a—y plane and sub- 
jected to the rf. field #,+i#,—EH,e' and the steady field H,, the 
currents induced by the rf. field will flow parallel to the surface, and 
the surface impedance Z, will be given by the expression 


Z2=4rtwalog is =e, tee ee 
(Gaussian units) so that 
Z, = (4arafogc?)4 ity 1) 
(b) If the applied r.f. field, instead of being circularly polarized, is 
arranged so that either H(0) or J, is linearly polarized, the associated 
surface impedances Z, and Z, can be obtained from Z, using equations 
similar to (2) and (3) respectively, with o, p replaced by 1/Z, Z; the 
resulting expressions agree with those given by Dresselhaus, Kip and 
Kittel (1955 b) and Donovan and Sondheimer (1953). (See also Donovan 
1954, 1955.) 


(c) If the specimen is arranged with its surface in the a—z plane, we can 
no longer assume the r.f. fields and currents to be circularly polarized 
in the a—-y plane, since we must now have J ,=0 (apart from the negligible 
currents needed to set up the space charges producing the Hall field £,) ; 
we must consider a linearly polarized current J, and apply eqns. (4). 
We then have Z?=47iwp,/c?, so that, from (4), Z is independent of 
wo and H, entirely. These results have already been obtained and 
discussed briefly by Anderson (1955). 

Now in attempting to observe cyclotron resonance experimentally 
conditions are normally chosen so that both w and 7 are as large as 
possible, since only if wz is large compared with unity can the resonance 
point wo=w be located accurately. But in any reasonably good metallic 
conductor the electronic mean free path will then be much greater than 
the skin depth, so that conditions are ‘ anomalous’ rather than classical, 
and eqns. (1)-(6) are no longer applicable. The problem must now be 
reformulated taking the finite free path into account. This is readily 
done for a quasi-free electron model using the ‘ kinetic’ approach 
(Chambers 1952) and it turns out that for the geometry of case (c) above 
the problem is quite insoluble analytically, but that cases (a) and (6) 
are more tractable ; for case (a) in particular the integral equation to be 
solved is formally identical with that occurring in the ordinary theory of 
the anomalous skin effect (Reuter and Sondheimer 1948, eqn. (18)), 
except that the quantity a=w7 occurring in that equation has to be 
replaced by w’7.¢ The variation of Z, with w’z is thus identical with the 
variation of Z with wz in ordinary anomalous skin effect theory, for which 


+ After completing this paper, it was found that Azbel’ and Kaganov (1954) 
had arrived at the same formulation ; these authors however discuss only the 
behaviour in the classical and extreme anomalous limits. 
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series expansions and tables have been given by Dingle (1953). His 
tables cover only positive values of w7, because only these are relevant 
for H=0, but in cyclotron resonance negative values of w’r must also be 
considered. His series expansions (6.3) and (6.4) have therefore been 
used to calculate values of Z, for w'r <0. 

The resultant curves are shown in the figure, where R,/R, and X,/X, 
are plotted against —w’r=(er/m*c)H,—wr. Here Z,=R,+iX,, and 
Ry, Xo are the values of R ,, X, for w'7=0. Curves are shown for several 
different values of the parameter « introduced by Reuter and Sondheimer : 
if d=(c?/27wo9)!” is the classically calculated skin depth, neglecting 

_ relaxation, and / is the electronic mean free path, «3/7/2562. Thus the 
classical limit corresponds to «<1 and the extreme anomalous limit to 
a>1. Now for most pure metals at 4°K and at microwave frequencies, 
~~104-108; it will be seen from the figure that under these conditions the 
variation of &,/Ry with —w’rz is very much less marked than in the classical 
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The variation with —w’r of R,/R, and X,/Xo, for different values of the 
parameter «. 


limit «—0, and extremely accurate measurements would be needed to locate 
the ‘resonance’ point wye=w at all closely. If the applied r.f. field is 
linearly rather than circularly polarized, corresponding to case (0) above, 
the values of Z, and Z, can again be obtained from Z, by using equations 
similar to (2) or (3); they vary even more slowly with w’r. These results 


462 R. G. Chambers on Cyclotron Resonance 


thus confirm Kittel’s (1954) suggestion that failure to observe cyclotron 
resonance in metals might be due to the large values of //5 encountered. 
If we write « in the alternative form 


vu 37 og (wr)? 
tig | Bez iid 


where v is the velocity of electrons at the Fermi surface, and put o//~10°° 


Gaussian units, v~108 cm/sec, representative of the values occurring 
in metals, we find «~1027(w7)?/w?. Now from the figure, we need w7=2 
and #< 100 say if experimental results are to lead to any useful estimate 
of m*: we thus require w> 10! sec!, or a wavelength AS200pu. And 
at w=10sec-1, if m*~m, the resonance condition wy—w requires a 
steady magnetic field H,~5 10° gauss. On this basis, therefore, there 
seems to be little hope of cyclotron resonance becoming a useful tool in 
the study of metals. 

We now have to consider how far these results are likely to apply to 
real metals. For this purpose it is helpful to introduce resistive and 
reactive skin depths 6, and 6,, such that in the classical region, where 
E(z)=E(0) e-*. 1/5,+7/5,=k. In the anomalous region, H(z) is no longer 
exponential, but we can define 5, and 6, by 


1/8, +4/8 = — (dB/dz))/E(0)=477iw|c2Z,, 


so that 8, and 5, can be derived from the measured surface impedance. — 


Now for the free-electron model in the ‘almost classical’ region, 
the departure of R, from its classical value is proportional to 1//6,., 
and that of X, to 1/6,\/(1+w”7?). It is plausible to assume that 
similar results will hold for more complex models or for real metals, and 
that R, or X, will begin to depart from classical behaviour when a 
suitably averaged value of / or l/\/(1+w’7?) becomes appreciable 
compared with 6, or 6;._ In the anomalous region, where these ratios are 
large, we may expect that, as for the free-electron model, R, and X, 
will vary with H, much more slowly than predicted classically. Thus, 
for example, minority carriers of small m* will probably produce a much 
less clear-cut perturbation of R&, than the classical theory (Anderson 1955) 
predicts. 

More complex models and real metals will differ from the free-electron 
model in three main respects. First, because of the spread of values of 
m* and the presence of holes as well as electrons, different parts of the 
Fermi surface will * resonate ’ at different values of H,, resulting in an even 
greater blurring-out of the curves. Secondly, in any but the simplest 
models v, for a given electron will vary periodically with time, so that the 
angle between the direction of motion and the surface of the metal will 
vary periodically. This makes the anomalous skin effect problem insoluble 
analytically, even for a simple spheroidal Fermi surface, but the 
‘ineffectiveness concept ’ approach (Pippard 1947, 1954) suggests that 
it may tend to increase both R, and X, somewhat when | wg7| becomes 
comparable with unity. Thirdly, any more realistic model will give a 


——— 
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non-zero d.c. magneto-resistance effect, and similarly the classically 
calculated values of Z,, and of 5, and 8;, will tend to increase with 
| 7.| more rapidly than for the free-electron model. For most metals, 
however, we shall still have 8,</ up to the highest values of lid | 
likely to be used, so that the metal will still be well in the anomalous 
region ; the tendency to approach the classical region as | H ,| increases 
will merely confuse the experimental results still further. For bismuth, 
however, the d.c. magneto-resistance is extremely large, and the results 
of Galt et al. (1955) suggest that a field | H,|~1 kg may increase 8, 
sufficiently to restore approximately classical conditions. Working 
at 4°K and 24 000 Me/s, they found that R, rose rapidly from its zero- 
field value as | H .| was increased up to about 1 kg, and thereafter varied 
much less rapidly and followed a course resembling that for a classical 
free-electron model. Now at 4°K and 24000 Me/s, the zero-field value 
of 1/5, for Bi is probably ~100-1000, corresponding to «~104—108 
(cf. Pippard and Chambers 1952, Chambers 1953); clearly we must 
suppose the effective value in fields >1 kg to be much smaller, to account 
for the observed results. Thus it is thought that the ‘trough’ in R, around 


zero field is a magneto-resistance effect; only for |H,|>1kg can the 


results be interpreted, tentatively, on a classical model, and then only 
if calculation shows that in this region //5, is not too large. Anderson 
(1955) has suggested an alternative interpretation of the trough, assuming 
the classical eqn. (5) to hold, in terms of a small contribution to o, from 
minority carriers of very small m*. It is thought that this interpretation 
must be rejected, however, not only because in zero-field conditions 
(5) is certainly not valid, but also because, even if (5) were valid, a peak 
or a ‘ dispersion curve ’ would be obtained rather than a trough, unless m* 
for the minority carriers exceeds m* for the majority carriers, which is 
quite incompatible with Anderson’s interpretation. 

It has been suggested (Dyson, private communication) that a more 
pronounced resonance might be obtained in metals under anomalous 
conditions if measurements could be made on a film of thickness ¢ small 
compared with the skin depth (and so, a fortiort, small compared with /), 
instead of on bulk metal. Jiven if such extremely thin films could be 
reliably produced, however, it is doubtful whether they would in fact 
show a very pronounced resonance : for the limiting case t<6</, where 
we may assume # constant across the film, it is easily found that for 
a free-electron model Z, is given by 


7 = e[m(z)+ouestac..| pe es. <77) 
K 


Z 4] 


where c=t(1--iw’7)/1 (ef. Fuchs 1938, Sondheimer 1950, eqns. (12), (25)), 
so that, for example, if #/J=10-*, R, increases by only 4% as |w’7 | 
increases from 0 to 1, and the minimum in R, at resonance is extremely 
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broad. We may note also that in the opposite limiting case of a classical 
thin film, wlth 1<t<8, we have simply 


Z = (lier fog 7 8) 


so that R, shows no resonance effects at all. 

Finally, it may be worth emphasizing the extent to which the effects 
observed in cyclotron resonance experiments will depend on the experi- 
mental arrangement. Any conductor placed in a static electric field 
normal to its surface will tend to screen out the field so that none appears 
in the interior. But if the conduction electron concentration is low enough, 
as in a semiconductor, the screening is imperfect, and if an r.f. field of 
sufficiently high frequency is used instead of a static field, it breaks down 
completely, so that the electric induction inside the conductor is equal to 
that outside. Thus if a small piece of the material is placed in a resonant 
cavity in a region of high electric field, the rate of power dissipation W 
will be proportional to #?R(c), and since the stored energy in the cavity is 
proportional to #?, the Q of the cavity (after correcting for other losses, 
including coupling losses) will be inversely proportional to R(c). If 
the concentration of carriers is higher, their effect in screening out the 
applied field can no longer be neglected ; it has been discussed in terms of 
the ‘ plasma frequency ’ by Dresselhaus, Kip and Kittel (1955 b). For 
a metal, of course, the screening-out of a normally applied field E is 
practically perfect, even at the highest microwave frequencies, and if the 
only way of inducing a current to flow in a metal was by putting a small 
isolated specimen in a region of high electric field, it would indeed be 
impossible to observe cyclotron resonance in metals, even under classical 
skin effect conditions (Kittel 1954). But in fact currents can be perfectly 
easily induced by making the specimen part of the wall of the cavity, 
in a region of high r.f. magnetic field : the current flow is then everywhere 
tangential to the surface and ‘plasma’ effects do not arise (cf. also 
Anderson 1955). We then have Wocl?R(Z), or WocH?R(Z), where H, 
is the magnetic field at the surface of the metal ; since the stored energy 
in the cavity is proportional to H,?, we now have the corrected Q of the 
cavity varying inversely with R(Z). 

The fact that a well-defined resonance peak is predicted theoretically 
for semiconductors but not for metals, even on classical skin effect theory, 
is due simply to the difference in experimental arrangements: semi- 
conductors are measured under essentially constant-voltage conditions, 
so that WocR(c), and metals under constant-current conditions, so that 
Woak(Z). Ifit were possible to measure a semiconductor under constant 
current conditions, with WocR(p), no resonance would be observed 
(cf. (1); ef. also (8)); conversely, if a metal obeying (6) could be 
measured under constant voltage conditions, with WocR(1/Z), a well- 
defined though asymmetrical peak would be obtained at w'’t=—1/,/3. 
However, under anomalous conditions the peak broadens out and the 
position of the maximum shifts (e.g. to w’t~—9 at «=1000), so that 


i) 
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even if such measuring conditions could be realized experimentally, 
cyclotron resonance experiments on metals would remain of doubtful 
value. 
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SUMMARY 


The yields of 7~ and z+ mesons from calcium have been measured at 90° 
to a photon beam with a maximum energy of 330+10 Mev. The mesons 
were detected in nuclear emulsions. The ratio of the numbers of 7 to 7* 
mesons was determined to be 0:86-0:18 at a pion energy of 80 Mev. 
Above 95 mev the ratio fell rapidly to 0:39--0-12 at 112 Mev. Approxi- 
mate values for the cross sections for the production of 7* and 7~ mesons, 
as a function of pion energy, were also determined. 

The experimental results may be given a semiquantitative explanation 
by the theory of Butler (1952). Preliminary calculations by Moorhouse 
and Laing (private communication) are in fairly good agreement with the 
measured data. 


§ 1. INTRODUCTION 


LirravER and WALKER (1952) have measured the yields of ~~ and z* 
mesons from a wide range of nuclei. Their results were obtained by 
observing pions with kinetic energies in the interval 50-80 Mev at an angle 
of 135° to a photon beam with a maximum energy of 310+10 Mev. The 
yield of negative to positive mesons (7 ~/7z* ratio) from calcium was- 
0-58-+0-06. This ratio showed the greatest deviation from unity for the 
symmetric nuclei investigated. Motz etal. (1955), using a linear accelerator 
to produce a bremsstrahlung beam with a maximum energy of 300 Mev, 
obtained a value of 0-70-L0-03 for the 7~/7+ ratio for calcium. In this 
experiment, mesons in the energy interval 43 to 63 Mev were detected at 
an angle of 75° to the direction of the photons. Both results are lower 
than the value calculated by Butler (1952) from a computation of the 
relative densities of neutrons and protons at the surface of the calcium 
nucleus. 

The present paper describes measurements made on the yields of 7~ and 
7* mesons from calcium as a function of their kinetic energy. Mesons 
emitted at 90° to the photon beam, with energies in the interval 60 to 
120 Mev, were detected in nuclear emulsions. The irradiating flux was a 
bremsstrahlung spectrum from the Glasgow electron synchrotron with 
photons of a maximum energy of 330-+10 Mev. 


+ Communicated by Professor P. I. Dee, F.R.S. 


Photoproduction of Charged Mesons from Calcium 467 


§ 2. EXPOSURE OF THE PLATES 

The plates used in the experiment were [ford G5 and C2 nuclear 
emulsions of 100 microns thickness and total area 4in.x2in. The 
geometry employed is shown in fig. 1; the top edge of each plate was 
6-4 cm from the centre of the calcium target and the copper absorber was 
2-20cem thick. To provide a homogeneous medium for stopping the 
mesons, eight plates were spaced at regular intervals in a slab of glass 
about 2 in. thick, and of the same composition as that used as a backing 
for the emulsions. Thus the emulsions were used to sample the number of 
mesons stopping in the glass as a function of distance from the target. 
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Experimental arrangement used to measure the yield of charged mesons 
from calcium at 90° to the direction of the photon beam. 


With this geometry, the top edges of the plates received mesons with a 
maximum angular spread of +37° and a probable angular spread of +11° 
about a direction at right angles to the photon beam, while for the bottom 
edges of the plates the maximum and the probable angular spreads were 
+15° and +5° respectively. The kinetic energies of the mesons were 
obtained from range-energy relationships for the materials traversed 
(i.e. calcium, copper and glass). Mesons in the energy interval 60 to 
120 mev could be detected with the apparatus. 


§3. COLLECTION OF DaTA 
The combined results of Menon et al. (1950) and Adelman (Cee) 
indicate that (48-7-0-9)% of all the 7 mesons ending in Ilford C2 
emulsions have two or more secondary tracks, and (23-3-0-7)% have 
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only one secondary track. Since all 7+ mesons ending in an emulsion 
produce single-pronged events (viz. 7—« decays), the ratio of 7 to mt 
mesons in a given number of events is simply 


oC 


n 
0-233 
0-487 Oj Gee 0-487 Cy 


where o,, is the number of mesons which end with two or more secondary 
tracks and o, is the number which end with only one secondary track. 
This method of separation overcomes the uncertainties involved in dis- 
tinguishing (a) . mesons from 7~ mesons with no secondary track, and (b) 
a— decays from 7~ mesons with only one secondary track. 

The plates were scanned for meson endings with one or more secondary 
tracks, and the position and nature of each event was noted, using the 
same convention for a star prong as that used by the authors referred 
to above. Several areas of each plate were scanned twice to determine 
the efficiency of scanning. The efficiency values obtained are indicated 
by comparing the figures given in columns 3 and 5 of the table with those 
in columns 4 and 6. 


(1) (2) (3) (4) (5) (6) (7) (8) 
62— 83 | 80-2] 146412 | 218426 | 6347-9] 70-54 9-3 | 0-784.0-16 | 0-86+0-18 
82— 89 | 87-6] 183414 | 246422 | 67+8-2| 76-6+10-1 | 0-75+0-14 | 0-84+0-16 
88— 98 | 94-2] 175413 | 237428 | 5747-6] 68-8+11-1 | 0-69+40-16 | 0-78+0-18 
97-106 | 103-0 | 104410 | 123415 | 2745-2] 31-8+ 6-5 | 0-60+0-16 | 0-69+0-19 
106-119 | 111-9} 91+10 | 109414 | 1443-7] 16-4+ 4-6 | 0-33+0-11 | 0-39+40-12 


(1) Energy interval (Mev); (2) mean energy (Mev); (3) observed number of one-prongs ; 
(4) corrected number of one-prongs (c,); (5) observed number of stars; (6) corrected 
number of stars (on); (7) uncorrected 7~/7* ratio; (8) a7/m* ratio corrected for absorption. 


The total yield of charged mesons at 90° (laboratory coordinates) was 
found as a function of meson energy, and is given in fig. 2. Due to 
uncertainties in the beam calibration, the absolute scale of the differential 
cross section has an estimated error of +30°%%. 

The data for the analysis of the 7~/z+ ratio are contained in the table, 
and the ratio as a function of pion energy is shown in fig. 3. At each 
point the errors given are statistical standard deviations, and include 
statistical errors arising from the efficiency determinations. The solid 
curves show the results when corrected for meson absorption. These 
corrections, although small (~10%), are somewhat uncertain. The 
method used is therefore given in an appendix. 


§ 4, Discussion 


The energy distribution of mesons obtained in the present work is 
consistent with a description in which the production arises from an 
interaction of the photon with a single nucleon in the nucleus. Butler 
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(1952) has shown that this production process, compared to that from a 
free nucleon, is modified as follows : 


(a) As a result of the spread in the momentum of nucleons in the target 
nucleus, there exists no unique photon energy corresponding to a given 
energy and angle of emission of the meson. Instead, for each energy and 


angle of emission of the meson there corresponds a band of photons avail- 
able for meson production. 


Fig. 2 
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‘Energy distribution of charged mesons produced from calcium at 90° by a 
“bremsstrahlung beam with a maximum energy of 330 +10 Mev. The 
solid curve shows the yield when corrected for absorption. 


(6) In general, for a given energy and angle of emission of the meson, 
the average energy of the photons available for production from a nucleus 
is different from the photon energy required for production from a free 
nucleon. The difference depends upon the binding energies of the target 
nucleons, and the interactions of the recoiling nucleons and mesons with 
the product nucleus. These interactions include the Coulomb interaction. 

The rapid change in the 7~/7+ ratio above 95 Mev can be explained in 
terms of the differing positions on the bremsstrahlung spectrum of the 
bands of photons available for the production of the 7~ and 7+ mesons 
respectively, and the existence of a cut off in the spectrum at 330 Mev. 
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As an illustration, we shall consider the factors influencing the pro duction 
of z- and z+ mesons emitted at 90° to the photon direction, with energies of 
90 and 110 Mev. In the case of meson production from the surface of the 
nucleus, Butler (1952) has estimated that the bands of photons are about 
50 Mev wide, and that they are centred (approximately) about the photon 
energy required for production from a free nucleon. We shall further 


Fig. 3 
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a/7* ratio as a function of pion energy for mesons produced from calcium 
at 90° by a bremsstrahlung beam with a maximum energy of 
330+10 Mev. The solid curve shows the ratio corrected for absorption. 


assume that, on the average, the photon energy required for the production 
of the 7~ meson is 10 Mev greater than that for the production of the z+. 
A difference in energy of this magnitude would appear reasonable when 
considering the differences in the binding energies of the last neutron and 
proton in the Ca nucleus. Thus az~ meson emitted at 90° to the photon 
direction with a kinetic energy of 90 Mev has available for its production a 
range of photon energies in the bremsstrahlung spectrum from about 
290 Mev to the cut-off at 330 Mev, while for a 7+ meson of the same kinetic 
energy and angle of emission the photon energies available are 280 to 
330 Mev. At a higher kinetic energy of the pion, say 110 Mev, the photon 
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energies available for the production of the ~~ and 7+ meson are respec- 
tively 320 to 330 Mev and 310 to 330 Mev. Thus a variation of the 7-/z+ 
ratio with pion energy of the type observed will result in this way. 

The value of the ratio at lower meson energies is in agreement with the 
value of 0-8 obtained by Butler (1952) by considering only the relative 
neutron and proton surface densities for calcium. The ratio obtained by 
Littauer and Walker (1952), as pointed out by Butler, is in a region where 
the effects of the cut-off in the photon spectrum cannot be neglected. 

Preliminary calculations by Moorhouse and Laing (private communi- 
cation), using an independent particle model with an oscillator potential, 
are in fairly good agreement with the experimental data for pion energies 
below about 100 Mev. The calculated results of these workers for the 
a /7* ratio are 0-79, 0-73 and 0-56 at pion energies of 80, 100 and 120 Mev 
respectively. Above about 100 Mev, the 7 /7z+ ratio becomes very 
sensitive to the shape of the bremsstrahlung spectrum and its maximum 
energy. 

It is hoped that an examination of the production at different angles to 
the photon direction will throw more light on the production process. If 
this process can be regarded as a modified free nucleon production, then 
the meson energy at which the total yield becomes small, and at which the 
7~/7* ratio rapidly varies, will be sensitive to the angle of production. 
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et Pek NED AL EX: 
The number of pions dN absorbed in a thickness dx of absorber is given 


by the relation dN 


ap no de ° 5 A . 5 ° A A (1) 


where n is the number of nuclei in unit volume of the absorber and c is the 
absorption cross section. The difference in the absorption of positive and 
negative pions is given by 


oN — EF = nfo ote ee hee), 
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where the superscripts refer to the sign of the meson. Bernardini ef ad. 
(1951), Bernardini and Levy (1951), and Ferretti and Manaresi (1955) have 
measured the mean free paths of 7~ and 7+ mesons in nuclear emulsions. 
The combined results of these workers indicate that the quantity («7 —o~) 
is independent of the kinetic energy of the meson in the range 10 to 90 Mev. 
Using this result, eqn. (2) may be integrated to obtain the relation 


No" ee exp £ = ae 
We Nee p {n(o-—ot)a} 
where NV, is the flux of mesons at the source. 

This relation was used to correct the observed 7 /7* ratio (i.e. the ratio 
N~-/N* in the above notation) to the ratio at the source. An average of 
the results of the authors given above yielded a value of 0-26 barns for 
(o-—o*t), where o is the cross section for star production together with 
inelastic scattermg by the nuclei of an emulsion, excluding hydrogen. 
Since the average geometrical cross section of the elements in the glass was 
0-67 of that for the complex nuclei in a nuclear emulsion, the quantity 
(o —o*) for glass was taken to be 0-67 x 0-26, to give 0-17 barns. For 
copper (o —ot)=0-46 barns, a value found experimentally by Cassels 
(reported by Burcham 1955). Both these values have large experimental 
errors (~30%). A change of 30°, however, in the value of (o~—a*) 
alters the value of the 7~/7* ratio by only 5%. 
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Iv is well known that there are not many examples of higher-order 
transitions which can be justifiably compared with Ehrenfest’s theory 
of such transitions. The specific heat C,, at a A-transition does not 
usually exhibit a well-defined discontinuity (the superconducting transition 
is exceptional in this respect), but rather rises on the low-temperature side 
to a very high maximum whose ultimate height appears to be limited only 
by the shortcomings of experimental technique. This behaviour, while 
diminishing the practical value of Ehrenfest’s method of approach, makes 
it possible to derive simple relations between various thermodynamic 
quantities which should be applicable in the vicinity of the \-transition. 
The entropy, considered as a function of temperature and pressure, 
rises steeply to a vertical or nearly vertical tangent along the A-line, and 
the analytical form of the entropy surface is dominated by this cireum- 
stance. In the neighbourhood of the A-line the surface may be regarded 
as cylindrical in the sense that the only principal curvature of significance 
is that in plane sections normal to the A-line. If then we assume that over 
short distances the A-line is straight and the entropy on the d-line is a 
linear function of temperature, we may write for S (p, 7’) the equation of a 
cylinder of arbitrary section whose generators are parallel to the A-line, 


Le: 8 (p, 1 )=S,+el+fi(el'—p), . . . «=. (1) 
in which S, and a are constants, f is a function describing the nature of the 


approach of the entropy to the A-point, and « is the pressure coefficient of 
the A-point, (dp/dT),. From (1) it follows that 


028 A Oe ‘4 Cd arte 
aq ed RVicEre seek and apie } 


hence, by use of Maxwell’s relations, over a short temperature interval 


a) -arn(@), ot (orn 
(5) G8), (Br (8) = 8 


in which f is written for 
1 (av 
Viet) 52 


the volume coefficient of expansion. From (2) we see that over such 
regions as the cylindrical approximation is justified (and it becomes 
oi ss URL OS ll te coe a 
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increasingly more exact as the A-line is approached) C,, and £ are linearly 
related, C,=aVT,B-+-const. > 0) eee) 


A similar argument applies to the volume, whose behaviour may be seen 
from (3) to follow that of the entropy closely, and yields 


B=aK-+ const. jee” & ee 


in which K is the isothermal compressibility. 


Fig. 1 
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It is easy to see that at a true second-order transition, where the slope 


of the entropy surface changes abruptly, (3) and (4) yield Ehrenfest’s 
relations for the discontinuities in C,, Band K, It should be mentioned 
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that Kuper’s (1955) analysis of the A-transition implies (3) and (4); in 
fact there is no particular novelty attached to these equations, except 
that they do not appear to be derived explicitly anywhere in the literature. 

The work of Simon, von Simson and Ruhemann (1927) and of Lawson 
(1940) on ammonium chloride provide suitable data for checking (3). 
In fig. 1 C,, is plotted against 8, values of the temperature being shown 
along the line (the transition temperature is 242-8°k). A remarkably good 
straight line is obtained, from which it may be inferred that the cylindrical 
approximation holds over at least four degrees below the transition 
temperature, that is, the shape of the specific heat curve is not highly 
pressure-dependent. From the slope of this line the temperature coefficient 


Fig. 2 
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T,v.p for NH,Cl. The line is drawn at the slope deduced from fig. 1. 


is deduced as 113 kg cm-2 deg-!. Figure 2 shows the transition tempera- 
ture at zero pressure and two transition temperatures at high pressures, 
according to Bridgman (1931), together with the calculated slope of the 
d-line at zero pressure. In view of the poverty of points the agreement 
may be regarded as satisfactory. Agreement with the data immediately 
above the transition temperature is not so good, the discrepancy being such 
as to suggest that the measured values of C,, are too high. This concords 
with an observation by Ziegler, quoted by Lawson (1940). The data 
for ammonium chloride have also been examined by Davies and Jones 
(1953) in connection with their thermodynamic theory of glasses, and been 
found not to be in accord with equations which they have derived on the 
basis of a model of the ordering process, and which bear a superficial 
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similarity to (3) and (4). In view of the good agreement found here with 
a purely phenomenological theory it seems likely that their discrepancy 
does not arise from the experiments so much as from the model or their 
interpretation of the parameters appearing therein. 

An attempt to check the equations by means of the A-point in liquid 
helium was not so successful. The values of C,, quoted in the literature 
are not accurate enough to be useful in this connection, so that only the 
validity of (4) can be examined. The data given by Atkins and Edwards 
(1955) on B and K appear good in themselves, but the value of « deduced 
from them differs by a factor of 2 from that measured directly. This 
discrepancy may be eliminated by assuming that the temperature scales 
for the two sets of data are shifted relative to one another by 0-001°K. An 
error of this magnitude is only too probable in practice. One may con- 
clude from this comparison that if the validity of (3) and (4) be accepted, 
the equations may provide a sensitive test of experimental data in the 
neighbourhood of a A-point. 
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ABSTRACT 


The saturation magnetic moments of certain iron alloys are difficult 
to interpret in terms of models found satisfactory for nickel alloys. It 
has been suggested that in pure iron the total number of empty d-band 
states is appreciably greater than the number of unpaired spins calculated 
from the saturation moment. The saturation moments of some iron— 
cobalt—aluminium alloys have therefore been measured and the effect 
of aluminium on the moments of iron—cobalt alloys is compared with its 
effect on the moment of pure iron. It is concluded that in iron the number 
of empty d-states exceeds the number of unpaired spins; the latter 
quantity is 2-22 per atom, while the value suggested for the former is 
greater than 3 per atom. The details of the effects found in the ternary 
alloys are discussed in terms of the d-band shape. 


§1. INTRODUCTION 


Ir is well known that when alloying elements such as copper, zinc, or 
aluminium are added to nickel the saturation magnetic moment at 
absolute zero falls as if atoms of the solute element contributed to the 
band structure of nickel a number of electrons equal to their conventional 
valency. Quantitative agreement with experiment results if it is assumed 
that the number of electrons per atom in the 4s-band (ns) is unaltered and 
that the electrons entering the empty states in the 3d-band all have spins 
opposed to those which give rise to the resultant moment in nickel. 
The latter assumption amounts to identifying the number of unpaired 
spins per atom (q) with the total number of d-band holes, that is, to 
assuming a condition of maximum multiplicity to obtain in the d-band 
(fig. 1(a)). 
When aluminium is added to iron, however, the effect on the saturation 
moment seems at first sight to be one merely of dilution, the moment per 
iron atom being unaffected. It has been suggested (Stoner 1947) that in 
iron the exchange interaction energy might be such that the separation of 
the half-bands containing spins parallel to and antiparallel to the resultant 
moment is insufficient to maintain maximum multiplicity (fig. 1 (6)). 


+ On leave of absence from Imperial College, London, 8.W.7. 
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Empty states would then exist in both half-bands, an alloying element 
such as aluminium could contribute electrons to both (with little resultant 
change in the moment per iron atom), and the number of unpaired spins 
per atom (q) given by the saturation moment at absolute zero would be 
less than the total number (nq) of d-band holes. 


Fig. 1 


(a) (b) 


Hypothetical d-band configurations (density of states plotted against energy) 
for ferromagnetic transition metals ; shading indicates occupied states. 
(a) Saturation moment g (in Bohr magnetons/atom)=total number of 
empty states per atom, ng. 
(6) ¢<na. 


The possibility of the ratio (€9) of q to nq being less than 1-0 has been 
discussed by Wohlfarth (1949 a) in connection with the nickel-copper 
alloys which have small interaction energies. This situation may also 
arise when the interaction energy is large if nq is large. Additions of 
cobalt to nickel increase nq but an examination of magnetic data for this 
system (Wohlfarth 1949 b) suggests that ¢, is unity over the whole compo- 
sition range. Saturation moment data for face-centred cubic cobalt- 
copper alloys (Crangle 1955) also indicate that ¢, is unity in this modifica- 
tion of cobalt. Additions of iron to cobalt still further increase the value 
of nq but an analysis of the type used for the nickel-cobalt alloys is rendered 
impossible by the change from body-centred cubic to face-centred cubic 
structures with increase of temperature. It has been suggested, however 
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(Hume—Rothery and Coles 1954) that the decrease of g with iron content 
in the iron-rich iron—cobalt alloys indicates a fall of €, from unity. It was 
therefore of interest to examine the effect of aluminium on the saturation 
moments of iron—cobalt alloys, and the results of such an investigation are 
presented below. | 


§ 2. MATERIALS, METHODS AND RESULTS 


Five iron—cobalt alloys were chosen as suitable materials to which 
aluminium additions might be made; these had cobalt contents of 
approximately 20, 35, 45, 50 and 60 atomic percent cobalt. Pairs of 
iron—cobalt—aluminium alloys were then prepared, each pair having an 
iron/cobalt ratio corresponding to one of the above compositions ; one 
of each pair contained about 5 atomic percent of aluminium and the other 
about 9%. The three pure metals used for the preparation of the alloys 
were of a purity greater than 99-9°,, and weighed quantities of them were 
melted together under argon in a high-frequency induction furnace. 
Homogenizing heat treatments of the ingots were carried out at about 
1350°c. The alloys were all too brittle for hot-forging to be possible, one 
of the ingots shattering in preliminary attempts. Small cylindrical 
specimens (about 8mm long and 4mm in diameter) were therefore 
machined from the ingots, and these were subjected to further homo- 
genization followed by a short annealing treatment in the body-centred 
cubic phase field above the ordering temperature of the binary iron—cobalt 
alloys. Microscopical examination of the specimens showed them to be 
single-phase after the above treatment. (The phase constitution of the 
ternary Fe-Co-Al system has been examined by Késter (1933) and 


pe 


Nominal Analytical Results Saturation Moments 

Atomic ar 
patio Atomic Atomic Bohr Magne- Eis ar ie nD 
Co/(Fe+Co) | Wt % Fe | Wt % Al . Ratio ; Fraction | tons per Atom | <0, metal 
Co/(Fe+Co) x of Al qd Ri oen g/(123) 

0 2-41 2-417 

0-198 { 77-0 3:05 0-198 0-0616 2-24 2:39 

L 75:9 4-70 0-197 0-0936 2°12 2-34 

ae 0 2-44 2-444 

oe el: 62-0 3-15 0-346 0-0638 2-23 2-38 

| 0 2-43 2-434 

0-441 53-0 2°85 0-442 0-0585 2°22 2-36 

| me Ogal 4-80 0-441 0-0961 2-09 2°31 

‘te ® 0 2-39 2-39T 

0-520 44: 2-70 0-525 0-0551 2-23 2-36 

: { 45-2 4-25 0-515 0-0864 2-08 2-28 

0 2-30 2-30T 

3 9.95 

0-610 36:8 2:35 0-610 6-0483 2-14 2-25 

{ 36-0 4-10 0-611 0-0837 2-00 2-18 


+ From the data of Weiss and Forrer (1929) 
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Edwards (1941).) Compositions :obtained by chemical analysis of the 
actual specimens are shown in the Table, together with the atomic 
fraction x of aluminium and the cobalt content expressed as the atomic 
fraction of the combined iron and cobalt contents ; this figure is the 
atomic fraction of cobalt in the binary iron—cobalt alloy to which the 
aluminium additions have effectively been made. 


Fig. 2 


QC Atomic fraction Al 
0:10 0-20 ©:30 


The effect of aluminium additions on the saturation moments of iron and 
binary iron—cobalt alloys ; x is the atomic fraction of aluminium, and 
q/(1—a) the saturation moment in Bohr magnetons per transition metal 
atom. The atomic percentage of cobalt in the binary alloy is indicated 
for each series of ternary alloys. The data for the iron-aluminium series 
are those of Fallot (1936) and Sucksmith (1933) 
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The saturation moment measurements were made at room temperature 
by a simple force method, the sample being suspended in that region of 
the pole-gap of an A. D. Little magnet where the field gradient is large 
and uniform over a vertical distance of about 0-5em. All measurements 
were made relative to a specimen of highly pure Armco iron for which the 
saturation moment was accurately known.+ Specimens could be 
exchanged without any displacement of the balance and suspension, so 
that there was negligible uncertainty in the reproducibility of the specimen 
position. Fields of up to 15 000 oersteds were employed and no variation 
of magnetization with field was detected for fields of more than 7000 
oersteds. All the alloys have Curie temperatures slightly higher than 
pure iron (Késter 1933) and only a very small error was introduced by 
assuming that the saturation intensities at room temperature bear the 
same ratio to those at absolute zero as does that of pure iron. If the same 
reduced magnetization against temperature curve applies to both iron 
and the alloys the ratio of the room temperature moment to that at 
absolute zero is 0-981 for pure iron and 0-987 for the alloys. From the 
measured saturation moments per gram and the known compositions 
values were derived for the saturation moment in Bohr magnetons per 
atom q and for the moment expressed as Bohr magnetons per transition 
metal atom g/(1—wz); these values are given in the table and values of 
q/(1—x) are plotted against the aluminium content in fig. 2. 


§3. Discussion 


In the earlier stages of the following discussion it is assumed that the 
general form of the 3d-band of the body-centred cubic iron—cobalt alloys 
is little affected by variations of the iron/cobalt ratio or by small additions . 
of aluminium, the main effect of compositional changes being to vary the 
number of occupied states in this band. The slight changes to be expected 
in the total d-band width and in the relative positions of the d-and 
3-bands do not affect the arguments of §§ 3.1 and 3.2 and are shown in 
§ 3.3 to be compatible with the model proposed. Discussion of the effects 
of variation of the exchange energy is also postponed to § 3.3. 


3.1. Saturation Moments and Total Numbers of d-Band Holes 


When copper is added to nickel one may regard the band system of the 
resulting alloy (overlapping 3d- and 4s-bands) as a property of the alloy 
as a whole, the effect of the extra electron possessed by each copper atom 
being to increase the number of occupied states. Since the density of 
states in the d-band is very much greater than that in the s-band nearly 
all the added electrons will occupy states in the former. When aluminium 
is added to nickel the behaviour of the saturation moment suggests that 
each aluminium atom contributes 3 electrons to the band structure, but 
it does not here seem justified to associate the d-states with the solute 


+ Its value was within 0-2°% of that given for iron by Weiss and Forrer 
(1929) in their study of the iron—cobalt alloys. 
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atoms. The most appropriate method of indicating the d-band configura- 
tion in such an alloy of a transition metal with a non-transition metal 
would therefore seem to be in terms of the number of unoccupied d-states 
per transition metal atom, and this quantity will be given (if one ignores 
changes in ns, the number of s-band electrons per atom) by the expression 


[na(1—2x) —a(v—ns)]/(1—2), 


where nq is the number of empty states per atom in the d-band of the pure 
transition metal, x the atomic fraction of the non-transition metal, and 
v its conventional valency. This expression will also be appropriate 
when additions of a non-transition metal are made to a binary alloy of 
two transition metals, nq and ns then referring to the binary alloy. 

As mentioned in §1, the above expression describes the behaviour of 
the saturation moments of many nickel alloys very well if it is assumed 
that qg, the saturation moment in Bohr magnetons per atom is equal to the 
total number of d-band holes. This assumption (that there are no paired 
holes) does not give a satisfactory description of the behaviour of iron 
alloys, as is clearly shown by fig. 2 where the curves marked N show the 
values given by it for g/(1—z), the saturation moment in Bohr magnetons 
per transition metal atom. Additions of aluminium to pure iron produce 
only a very small change in q/(1—2), and if the tentative interpretation 
of §1 (that involving a band configuration of the type shown in fig. 1(d) ) 
is correct, electrons contributed by the aluminium atoms would seem to 
be entermg the two half-bands in almost equal numbers. At larger 
aluminium contents, however, marked changes should be observed when 
all holes in one half of the d-band have been filled. 

In fig. 2 the curves marked P show the behaviour to be expected of 
q/(1—a) if the values of nq for the iron—cobalt alloys were those shown in 
fig. 3, and if electrons contributed by the aluminium entered the two 
halves of the d-band in equal numbers. (The latter assumption would 
require the densities of states in the two half-bands to be the same.) 
The experimental results do not show an initial constancy of qg/(1—2) 
for all cobalt contents, and the densities of states in the two halves of the 
d-band evidently differ in the alloys. This might have been expected, 
for on the collective electron treatment a rectangular band is inconsistent 
with 0<q/na<1-0 (Wohlfarth 1951). The final rate of fall of g/(1—2) 
with increasing # is, however, independent of the band shape, since it is 
characteristic of aluminium contents sufficient to fill all holes in one half- 
band, and for the alloys where this region is reached good agreement 
between theory and experiment can evidently be obtained by using the 
nq Values given in fig. 3. The corresponding values of ns and the observed 
values of g are also shown in this figure. The nq values were chosen to give 
agreement with experiment for the effects of aluminium on the 52% 
and 61%, cobalt alloys and to pass through the value of g obtained for pure 
cobalt by extrapolation from the body-centred cubic alloys. This value 
(1:80) is only slightly greater than that (1-75) for face-centred cubic 
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cobalt (Myers and Sucksmith 1951, Crangle 1955), and evidence that q is 
equal to mq for that modification of cobalt has been referred to in Sale 
The iron-rich portion of the ng curve is quite uncertain, but a lower limit 
of 2-9 can be set for the nq value of pure iron from the aluminium content 
of iron—aluminium alloys at which ¢/1—x) begins to fall sharply (Fallot 
1936, Sucksmith 1939). This is, however, only a lower limit, for this 
fall takes place in the composition range where ordering effects become 
apparent. | 
Fig. 3 


3°2 


Dee 


Je Nee 2-0 


20 40 6 
Fe, - Atomic % Co 


Saturation moments in Bohr magnetons/atom qg and postulated numbers 
per atom nq and ng of d-holes and s-electrons in the body-centred iron— 
cobalt alloys. 


S 80 Co 


3.2. Band Shapes and Alloying Effects 


Some indication of the relative magnitudes of the densities of states 
N(Z) in the half-bands (p and a) containing electrons with spins parallel 
and antiparallel to the resultant moment can be obtained from the 
experimental results for the alloys with 20%, 35%, and 44%, cobalt, 
for these relative magnitudes will govern the g/(1—) values of alloys with 
aluminium contents insufficient to bring about filling of the p half-band. 
For g/(l—x) to be constant in such alloys Np(#) must equal N,(#) ; 
if it is smaller g/(1—z) will fall with increasing 2, and if it is larger q/(1—2) 
will increase. The results for the 20%, 35%, and 44% cobalt alloys 
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indicate that N,,(e) is less than NV,(#) ; and this is not surprising, for the 
p half-band is more nearly full. For the 52% and 61% cobalt alloys the 
experimental results seem to correspond to aluminium contents sufficient 
to bring about filling of the p half-band. The simplest band shape which 
would account for this behaviour is that shown in fig. 4(a), where the 
Fermi level in the half-bands of pure iron is indicated by F,, and that in 
the iron-rich iron—cobalt alloy by F,. For convenience a single diagram 
is used here and in fig. 4(b) to represent the configurations at two different 


Fig. 4 


(2) (») 


a half-band 


p half-band 


Hypothetical forms for the d-band density of states curve in ferromagnetic 
alloys. #, is the Fermi level in an alloy where N,(£)=N,(Z) and F, 
that in an alloy where N,(#)>N,,(2). 


compositions. It is not intended to imply that the exchange energy is 
composition-independent. In view of the indications of a double- 
peaked d-band structure given by theoretical calculations (Fletcher and 
Wohlfarth 1951, Slater and Koster 1954, Koster 1955) consideration 
could also be given to configurations of the type shown in fig. 4 (b). 
An interesting consequence of the latter configuration is that additions 
to iron of manganese, which would increase nq and lower the Fermi 
level, might produce alloys in which q/(1—.) would rise for initial additions 
of aluminium. 
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3.3. Interaction Energy and Band Widths 


As pointed out by Wohlfarth (1949) in his discussion of the nickel— 
copper alloys the ratio (C9) of q to nq can fall below unity when the exchange 
interaction energy becomes small. The parameter which governs ¢, 
in the collective electron theory of ferromagnetism (Stoner 1947) is the 
ratio of interaction energy to the energy range of unoccupied states in the 
d-band. Consideration of the latter quantity shows that a fall of Gr 
as mq increases. in the iron—cobalt alloys need not be ascribed to any 
marked fall in the interaction energy. On the basis of a number of 
general arguments considerable increases in the energy range (AZ) 
of unoccupied states are to be expected as the iron-content of iron—cobalt 
alloys is increased. 


(a) If no change in the relative positions and breadths of the 3d- and 
4s-bands took place between nickel and cobalt and between cobalt and 
iron mq would increase by a little less than 1-0 in each step (a fraction of 
the removed electron coming from the 4s-band) and 4H would show a 
marked increase. 


(6) On general grounds it is to be expected that the overall width of. 
the 3d-band will increase through the series Cu-Ni+Co+Fe+Mn-Cr, 
and for a given number of d-holes 4 can be expected to be larger in iron 
than in nickel. 


(c) Atomic 3d levels move to higher energy relative to 4s ones along 
the above series of elements, and corresponding relative movements are 
to be expected in the metallic state. Such movements would lead to 
increases in mq of more than 1-0 holes per atom in the steps from nickel 
to cobalt and from cobalt to iron, and more rapid increases in 4# than 
suggested by (a) ; .and the nq values of 0-6 for nickel and 1-75 for face- 
centred cobalt show the expected trend. (The postulated values of 1-8 
and ~3-2 for body-centred cubic cobalt and iron are in accord with such 
movements.) 


(d) If the double-peaked structure for the d-band shown by calculations 
is correct, increases in nq will produce larger increases in 4H when nq 
is ~2-0 than when vq is ~0-6, for the Fermi level for ng=2-0 is likely to 
lie at a region of the a half-band where the density of states is smaller than 
for na=0-6. It is significant that the electronic specific heat coefficient 
of iron is somewhat less than that for nickel, in spite of the contribution 
made in the former metal (if the present model is correct) by the p 
half-band. 


The above arguments have been brought forward to show that, for a 
given exchange energy, ¢, values of less than unity become more probable 
as the iron content incréases. Such values are still more probable if the 
exchange energy decreases with an increase of the iron content, and the 
behaviour of the Curie temperature is in accord with a variation of that 


type, 


486 On Saturation Moments and d-Band Configurations in Iron 


ACKNOWLEDGMENTS 


The authors’ thanks are due to many colleagues whose assistance 
contributed to the work described above; in particular to Dr. J. E. 
Goldman and Dr. A. Arrott for much stimulating discussion, and to 
Dr. M. Simnad and the staff of the Metals Research Laboratory for placing 
metallurgical facilities at their disposal. Grateful acknowledgment for 
financial support is made to the U.S. International Cooperation Adminis- 
tration, the U.S. Office of Naval Research, and the U.S. Army Signal Corps. 


REFERENCES 


CRANGLE, J., 1955, Phil. Mag., 46, 499. 

Epwarps, O.S., 1941, J. Inst. Met., 67, 67. 

Fatuot, M., 1936, Ann. Phys., Paris, 6, 305. 

FLEetTcHer, G. C., and WonirartH, E. P., 1951, Phil. M. ag., 42, 106. 

Hume-Rotnery, W., and Coss, B. R., 1954, Advances in Pha ysics, 3, 149. 

Koster, G. F., 1955, Phys. Rev., 98, 901. 

Kostir, W., 1933, Arch. Bisenhuttenwesen, 7, 263. 

MYERS, ae P., and SucxsmituH, W., 1951, Proc. Roy. Soc. A, 207, 427. 

SLATER, RAG ss and Koster, G. F., 1954, Phys. Rev., 94, 1498. 

STONER, E. GC, 1947, Rep. Progr. Phys., 11, 43. 

Sucxsmita, W., 1939, Proc. Roy. Soc. A, 171, 525. 

Weiss, P., and Forrer, R., 1929, Ann. Phys., Paris, 12, 279. 

WoutrartH, E. P., 1949a, Proc. Roy. Soc. A, 195, 434; 1949b, Phil. Mag., 
40, 1095; 1951, Ibid., 42, 374. 


Pads] 


L. REVIEWS OF BOOKS 


Vistas in Astronomy. Volume I. Edited by A. Brrr. (London : Pergamon 
Press.) [Pp. xvi+776.] Price £9 9s. 
Taz two volumes of Vistas are dedicated to F. J. M. Stratton in honour of his 
seventieth birthday. The exertions of the Editor and originator of the scheme 
have brought together contributions from over two hundred astronomers, 
geophysicists, physicists and mathematicians to produce “a generous and 
thorough cross-section through the whole of contemporary astronomy and 
allied sciences”. The first volume ranges from matters concerning inter- 
national organization in astronomy to radio astronomy and solar physics. The 
contributions vary from autobiographical records, to announcements of new 
discoveries and reviews of particular branches of astronomy. The result is a 
book possessing a wholly individual character of remarkable vitality. It will be 
read and treasured long after more orthodox handbooks have served their day, 
so keeping alive the memory of one whose influence has so widely pervaded the 
work of his own times. The Editor and Publishers deserve admiring congratu- 
lations for a great undertaking carried through to complete success. 
W. H. McC. 


Numerical Analysis with emphasis on the application of numerical techniques to 
problems of infinitesimal calculus in single variable. By Z. Kopau. (Chapman 
and Hall.) [Pp. xiv-+556.] Price 63s. 

THE scope of this book is described more exactly by the full title than by the 
short title “‘ Numerical Analysis ’’, since such topics as the solution of problems 
in linear algebra and the numerical solution of partial differential equations are 
not included. The book is a written version of a course of lectures given to 
students at the Massachusetts Institute of Technology, and the author has 
preserved his rhetorical questions and other lecture room mannerisms. This 
accounts for the great and, indeed, somewhat excessive length of the book. 
It also accounts, however, for the book’s most distinguishing virtue, namely 
that the material is presented with a lively and stimulating commentary. 
Thus, although the treatment of particular formulae and methods is less concise 
than in ordinary books on numerical analysis, it is also less superficial. Any 
practitioner in numerical analysis, on dipping into the book, is likely to find 
something which will set him thinking, whether it be to extend, to amplify, or 
perhaps in some cases to refute what Dr. Kopal has said. Of particular value 
are the bibliographical notes at the ends of the chapters which contain critical 
reviews of the literature. Coming at a time when the development of automatic 
digital computers has called into question much of what was fast becoming 
established doctrine among numerical analysists trained in desk computing 
techniques, this book will form a valuable addition to the available literature. 


MrV.oW. 


Introduction to Modern Physics. 5th Edition. By F. K. Ricurmysr, E. H. 
Kennarp and T. Lauritsen. (New York: McGraw-Hill.) [Pp. xv+ 666. ] 
Price 56s. 6d. 

A irra edition of the well-known text book by Richtmyer and Kennard is very 
welcome. Professor Richtmyer died in 1939 and Professor Lauritsen has 
joined Professor Kennard in preparing this revised edition. The book has been 
brought up to date by including an account of the newly discovered phenomena 
of high energy physics, while the chapters on X-rays and physics of solids have 
been somewhat cut down, 
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Solid State Physics. Advances in Research and Applications. Volume I. 
Edited by FrepErick Serrz and Davip TurNBULL. [Pp. xii+469.] (New 
York : Academic Press Inc.) Price $10.00. 


Iris planned to publish two of these review volumes each year until the subject 
is adequately surveyed and then to publish annually. The present volume, 
the first of the series, contains the following articles. 

Methods of the One-Electron Theory of Solids by J. R. Reitz; Qualitative 
Analysis of the Cohesion in Metals by E. P. Wigner; The Quantum Defect 
Method by F. S. Ham; The Theory of Order—Disorder Transitions in Alloys 
by T. Muto and Y. Takagi; Valence Semiconductors, Germanium and Silicon 
by H. Y. Fan; and Electron Interaction in Metals by D. Pines. 

“There is no doubt that the volume will be very useful indeed, covering as 
it does the important recent advances in the electron theory of solids and 
much of the work both experimental and theoretical on the semi-conducting 
properties of silicon and germanium. It is intended that later volumes will 
pay more attention to the experimental and technical aspects of the subject 
than this one does. 

It remains to ask, what is ‘Solid State Physics’ ? Not, one gathers from 
the Introduction, crystallography or metallurgy, yet it includes the study of 
defects in crystal structures and the theoretical study of the binding forces 
in metals. Certainly a large part of the subject is concerned with conduction 
in solids. Perhaps it is best to define the subject by naming some of the people 
who work in it and who so describe their work, physicists who investigate the 
properties of solids without wishing to bind themselves to the narrower 
disciplines of metallurgy and crystallography. ' N. F.M. 


BOOK NOTICES 


Topological Transformation Groups. By D. Montgomery and L. Zrppi. 
(New York and London: Interscience Publishers.) [Pp. xii+282.] Price 
$5.50. 


The Number-System. By H. A. Taurston. (London and Glasgow: Blackie 
& Son Limited.) [Pp. viii+134.] Price 30s. 


Hydrodynamic Stability. By C. C. Lix. (Cambridge: University Press.) 
[Pp. xi+155.] Price 22s. 6d. 

La Teoria Fisica. By A. M. Dell’Oro, (Padua: Cedam.) [Pp. viii+-155.] 
Price not stated. 


Ludwig Boltzmann. By Engelbert Broda. (Vienna: Franz Deuticke.) 
[Pp. vili+-152.] Price (cardboard bound) DM. 9.50, (clothbound) DM. 11. 


[T'he Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Event SQ 1437. The two tracks marked A and B apparently diverge from a 
common origin in the lead above the cloud chamber. If this interpre- 
tation is correct, both tracks are due to negatively charged particles 


of mass greater than 700m,. The measurements made on the tracks 


are given in the table. 
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Event SF 992. The two tracks indicated by arrows intersect in the copper 
plate across the middle of the cloud chamber. The near-vertical, 
dense track has a curvature corresponding to a momentum of 
204-+14 Mev/e and its ionization density is estimated to be between 
51) and 8/5. From these figures the mass of the particle responsible 
is between 920m, and 1450 m,. The event is interpreted as an upward- 
moving, positive K-meson which is brought to rest in the plate and 
decays there. The secondary particle has a momentum of 154-20 Mev/c. 


